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Abstract 


SOLVING  A  MULTIACTIVITY  MULTIFACILITY 
CAPACITY-CONSTRAINED  0-1  ASSIGNMENT  PROBLEM 


by 

Krishan  Lai  Chhabra 

Richard  Martin  Soland,  Director  of  Research 


A  branch-and-bound  solution  algorithm  and  a  computer  program 
implementing  this  algorithm  are  developed  to  solve  a  multiactivity 
multifacility  capacity  constrained  0-1  assignment  problem.  The  math¬ 
ematical  formulation  for  such  a  problem,  called  problem  (P),  is  to  find 
x..  and  y,  values  that: 

ij 

m  n  p 

Minimize  Z  I  a . .  x . .  +  E  b.  y,  (i) 

i-i  j-i  13  13  k=i  k  k 


m 


subject  to  E  x  =  1 

i=l  J 

j =2 ,  •  •  • >  n 

(ii) 

m  n 

l  Z  d  x  < 

i-i  j-i  ijk  13  “ 

s.  Yu 
k  k 

k=l, . . . >p 

(iii) 

x .  ,  =  0  or 
ij 

1  for 

all  i  and  j 

(iv) 

yk  =  0  or 

1  for 

all  k 

(v) 

ii 


where  i  ,  j  >  k  are  indices  for  designs,  activities,  and  facilities, 
respectively;  x_  has  value  1  if  and  only  if  activity  j  uses  design 

i  ,  and  y^  has  value  1  if  and  only  if  facility  k  is  used.  A  design 

involves  the  use  of  one  or  more  facilities,  and  the  same  design  may  be  used 
by  several  activities. 


Problem  (P)  has  the  objective  of  minimizing  the  sum  of  a.^s  — 

the  variable  costs  due  to  the  assignments  of  activities  to  designs,  and 
b^ 1  s  —  the  fixed  costs  due  to  the  facilities  used.  Constraints  (ii) 


and  (iv)  ensure  that  each  activity  is  assigned  to  a  single  design. 


ijk 


is  the  capacity  required  at  facility  k  if  activity  j  uses 


Each 

design 


i  ,  and  is  thus  equal  to  zero  if  design  i  does  not  involve  the  use  of 
facility  k  .  Constraints  (iii),  therefore,  ensure  that  for  each  facility 
k  used,  the  total  capacity  required  does  not  exceed  the  capacity  available 
,  The  difficulty  in  solving  problem  (P)  stems  from  the  indirect  relation¬ 


ship  between  the  assignments  and  facilities, 
bears  on  all  the  constraints  (iii)  for  which 
therefore,  on  several  y^  variables. 


l .  e  . , 


ijk 


an  assignment 
is  positive. 


and , 


1 


The  branch-and-bound  solution  algorithm  uses  Lagrangian  relaxation 
as  a  basic  step  in  obtaining  lower  bounds.  In  addition,  it  includes  several 
operational  rules,  such  as  a  branching  rule  for  a  judicious  choice  of  the 
branching  variable,  a  capacity  rule  to  eliminate  infeasible  assignments, 
and  a  bounding  rule  to  eliminate  non-opt imal  assignments. 


This  dissertation  includes  relevant  background  leading  to  the 
formulation  of  problem  (P) ,  mathematical  development  of  the  branch-and- 
bound  solution  algorithm,  a  detailed  test  example,  and  computational  test 
results  using  the  computer  program.  The  areas  of  application  are  identified, 
and  suggestions  for  further  improvement  of  the  branch-and-bound  solution 
algorithm  are  included. 


The  computer  program  has  been  written  in  FORTRAN  IV.  A  detailed 
description  of  the  computer  program  and  guidelines  for  its  use  are  included 
in  a  separate  document  entitled  "Program  Description  and  User's  Guide  for 
ZIPCAP — a  Zero-one  Integer  Program  to  solve  multiactivity  multifacility 
Capacity-constrained  Assignment  Problems. M  Although  developed  for  capac¬ 
itated  problems,  the  computer  program  can  also  be  used  to  solve  uncapacitated 
problems  in  which  it  is  assumed  that  the  facilities  have  infinite  capacity. 
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1.  INTRODUCTION 


Mul  tiac t  i vity  multifacil  ity  assignment  p?*obleins  arise  in  such 
diverse  areas  as  public  health  care  systems  and  private  multi-echelon 
inventory/distribution  systems.  Such  systems  involve  the  assignment  of 
activities  or  tasks  to  groups  of  facilities  in  such  a  way  that  total 
system  cost  is  minimized.  The  total  system  cost  has  components  (fixed 
costs)  that  depend  on  the  facilities  actually  used  as  well  as  components 
(variable  costs)  that  depend  solely  on  the  assignment  made.  Most  recently 
[Gross,  Pinkus,  and  Soland  (1979)]  there  has  been  interest  in  including 
facility  capacity  constraints  as  well.  For  this  kind  of  problem,  i.e., 
a  multiactivity  multifacility  capacity-constrained  0-1  assignment  problem, 
we  have  developed  a  solution  algorithm  of  the  branch-and-bound  type  and 
a  computer  program  based  upon  it. 

The  computer  program  and  guidelines  for  its  use  are  described  in 
a  separate  document  [Chhabra  and  Soland  (1980)]  titled  "Program  Descrip¬ 
tion  and  User's  Guide  for  ZIPCAP  —  a  Zero-one  Integer  Program  to  solve 
multiactivity  multifacility  Capacity-constrained  Assignment  Problems." 

This  document  describes  the  development  of  the  solution  algorithm 
and  compu ta t ional  test  results  using  the  computer  program.  Suggestions 
for  further  improvement  in  the  solution  algorithm  are  also  included. 

This  chapter  reviews  the  relevant  literature,  provides  background 
leading  to  the  mathematical  formulation  of  the  mul t iac t ivity  multifacility 
capacity-constrained  0-1  assignment  problem,  called  problem  (P) ,  and 
includes  potential  areas  of  application.  The  theoretical  base  for  develop¬ 
ing  the  algorithm/methodology  are  described  in  Chapter  2.  Various  components 
of  the  methodology  are  covered  in  detail  in  Chapter  3.  Chapter  4  provides 
an  overview  of  the  computational  procedure  and  the  computer  program,  whereas 
computational  test  results  are  given  in  Chapter  5.  Suggestions  for  further 
research  and  potential  improvements  in  the  algorithm  are  included  in  Chapter 
6. 


-  i  - 
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It  may  be  noted  that  the  basic  terminology,  described  below, 
in  the  formulation  of  problem  (P)  includes:  activities  that  must  be 
assigned,  facilities  which  serve  the  activities,  designs  involving  one 
or  more  facilities,  fixed  costs  associated  with  tne  facilities,  and 
variable  costs  associated  with  the  assignment  of  activities  to  designs. 

The  following  review  of  the  relevant  literature  starts  with  the 
classical  assignment  problem  and  leads  to  the  fonnulation  of  problem  (P) . 
Different  authors  have  used  various  terminologies  in  describing  relevant 
formulations.  In  the  following  discussion,  the  original  terminologies 
are  used,  and  are  followed  by  our  equivalent  terminology,  where  appro¬ 
priate,  shown  in  parenthesis. 

1 . 1  Generalized  Assignment  Problem 

In  a  classical  assignment  problem  [Hillier  and  Lieberman  (1980)], 
the  purpose  is  to  find  optimal  pairs  of  agents  ar d  tasks  or  activities. 
Each  task  is  assigned  to  a  single  agent,  and  each  agent  is  given  a 
single  task,  and  the  suitability  of  a  particular  set  of  assignments 
is  determined  by  a  single  criterion  function  such  as  minimization  of 
cost.  In  a  generalized  assignment  problem  (GAP),  several  tasks  can  be 
assigned  each  agent,  subject  to  the  resources  available  to  the  various 
agents  [Ross  and  Soland  (1975)],  e.g.,  assigning  software  development 
tasks  to  programmers  and  assigning  jobs  to  computers  in  a  computer 
network. 

A  variety  of  well-known  facility  location  and  location-alloca¬ 
tion  problems  have  been  shown  to  be  equivalent  to,  and  therefore 
solvable  as  GAP T s  [Ross  and  Soland  (1977)].  Here,  in  general,  the 
tasks  represent  demand  centers  for  a  good  or  service,  and  the  agents 
represent  supply  centers  to  be  established  at  potential  sites  or 
locations.  Each  demand  center  must  be  supplied  from  a  supply  center. 

A  fixed  cost  is  incurred  for  each  supply  center  established,  and,  in 
addition,  there  is  a  cost  incurred  for  each  unit  processed  at  a  supply 
center  and  transportation  costs  incurred  for  the  units  sent  from  supply 
centers  to  demand  centers.  The  problem  may  be  "uncopac itated"  —  when 
there  is  no  limit  to  the  number  of  units  that  may  be  processed  by 


a  supply  i'ouut,  or  "capacitated"  —  when  there  are  restrictions  on  the 
number  o 1  units  that  may  be  processed.  The  objective  is  to  select  supply 
center  locations  and  set  up  a  distribution  assignment  so  that  the  total 
cost  is  minimized. 

1 . 2  Mult iac tivity  Mul t i f ac il ity  Uncapacitated 
Assignment  Problem 

A  salient  feature  of  the  above  facility  location  problems  is  that 
eacli  demand  center  (activity)  is  assigned  to  a  single  supply  center 
(facility).  Sometimes,  however,  it  may  be  desirable  to  assign  an  activity 
to  more  than  one  facility.  This  leads  to  the  concept  of  "design,"  and 
the  multiactivity  multifacility  assignment  problem  [Pinkus,  Gross,  and 
Soland  (197J) J .  Before  describing  such  a  problem,  some  terminology  is 
considered  first. 

A  design  involves  the  use  of  one  or  more  facilities,  and  represents 
a  meaningful  configuration  of  facilities  along  with  a  meaningful  strategy 
for  using  them  —  as  illustrated  in  the  following  examples. 

Consider  five  facilities  and  their  locations  as  shown  in  Figure 
1(a).  (From  practical  considerat ions ,  these  may  be  existing  and/or 
potential  locations.)  Three  of  the  possible  designs  are  shown  in 
Figures  1(b)  to  1(d) ,  Design  1  is  completely  centralized  since  it  uses 
only  one  facility,  whereas  design  3  is  completely  decentralized  since  it 
uses  all  the  facilities. 


© 

© 

© 

©  i 

(a) 

Locations  for  five 
facilities 

© 

© 

© 

(a)  Design  2:  three  facili¬ 
ties — locations  ©.  ©, 

© 


© 

(b) 

Design  1:  one 

facility- 

locat  ion(3) 

© 

© 

© 

© 

©. 

(d)  Design  3:  five  facilities 
all  locat ions 


Figure  1.  Examples  of  alternate  designs 

for  a  system  of  five  facilities 
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It  is  possible  for  several  designs  to  have  the  same  facilities 
but  different  configuration  and  strategies  for  using  these  facilities, 
e.g.,  a  multiproduct  mui ti-echelon  inventory  system  [Gross,  Pinkus, 
and  Soland  (1979)].  Figure  2(a)  shows  design  1  containing  certain 
facilities  (warehouses)  at  the  central,  regional,  and  local  levels 
or  echelons.  Figure  2(b)  shows  design  2  with  the  same  facilities 
but  having  u  different  configuration. 

Level  cr  Echelon 
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(a)  Design  1  (b)  Design  2 


Central 

Regional 

Local 


F igur e  2 . 


Example  of  alternative  designs  having 
che  sane  facilities  but  different 
c  onf 1 c u r a t io  n 


The  distribution  of  a  given  activity  at  various  facilities  under  design  1 
would  be  different  than  under  design  2,  depending,  of  course,  on  the 
inventory  policies.  This  results  in  different  variable  costs  (described 
later)  for  chat  activity  under  design  1  as  against  design  2.  In  fact, 
it  is  possible  to  have  a  situation  where  two  or  more  designs  have  the 
same  facilities  and  the  same  configuration  but  different  strategies, 
resulting  in  different  variable  costs.  For  example,  one  strategy  might 
specify  an  equal  distribution  of  a  specific  activity  over  the  various 
facilities,  whereas  another  strategy  could  impose  a  different  distribu¬ 
tion  scheme  over  the  same  facilities. 

In  general,  if  a  system  is  to  be  composed  of  at  most  p 

facilities,  the  number  of  alternative  designs  is  2^-1  if  no  two  designs 
have  the  same  facilities.  However,  with  the  same  facilities  but  different 
configurations  and  strategies,  the  number  of  alternative  designs  could  be 
much  higher.  In  practice,  it  is  possible  to  eliminate  a  majority  of 
alternative  designs  because  of  geographical,  political,  economical,  and 
other  factors. 
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The  anil  t  iaotiv  ity  multifacility  assignment  problem  seeks 
minimi /.at  iuu  of  some  measure  of  total  system  cost  such  as,  total 
expected  cost  over  a  given  time  period  or  total  discounted  cost  over 
the  lifetime  of  the  system.  The  system  cost  will  Include  investment 
costs  for  building  or  leasing  the  system,  operating  costs  for  operation 
and  maintenance  of  the  system,  and  the  costs  for  providing  necessary 
services.  Both  the  investment  costs  and  the  operating  costs  have 
fixed  as  well  as  variable  components  [Ross  and  Soland  (1980)].  The 
fixed  components  include  those  costs  associated  with  the  facilities 
of  a  given  design  which  are  independent  of  the  activities  served. 

Such  costs  are  called  fixed  costs.  On  the  other  hand,  the  variable 
components  and  the  service  costs  include  those  coses  which  are 
completely  dependent  on  the  service  demand  of  the  activities  at 
the  various  facilities  in  a  given  design.  Such  costs  are  called 
variable  costs.  By  definition,  both  the  fixed  costs  and  the  variable 
costs  are  relative  terms. 

An  equivalent  formulation  of  the  multiact ivity  multifacility 
assignment  problem  defined  by  Pinkus,  Gross,  anc  Soland  (1973)  is  as 
follows. 


Let  a , 


ik 


variable  cost  of  activity  j  using  design  i 
( l-l , . . . , m,  j-1 , . . . , n) 

fixed  cost  of  facility  k  (k=l,...,p) 

1  if  facility  k  is  included  in  design  i  , 

0  otherwise. 
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The  objective  function  of  this  problem  consists  of  two  distinct 
parts.  The  first  part  represents  the  total  variable  cost,  and  the 
second,  the  total  fixed  cost  of  the  system.  Constraints  (2)  and  (3) 
ensure  that  each  activity  is  assigned  to  a  single  design.  Of  course, 
the  optimal  solution  may  involve  the  use  of  more  chan  one  design. 

Problem  (PU)  is  a  0-1  nonlinear  programming  problem  (because  of 
the  step  function  u  ),  and  a  branch-and-bound  algorithm  using  linear 
underestimates  for  the  nonlinear  part  of  the  objective  function  has  been 
described  in  Pinkus,  Gross,  and  Soland  (1973).  A  heuristic  procedure 
for  this  problem  is  given  by  Khumawala  and  Stinson  (1980)  in  an  unpublished 
paper.  This  procedure  is  an  extension  of  some  earlier  work  [Khumawala 
(1973)]. 

1 . 3  Adding  Cap  icity  Constraints  —  Pro blem  (P) 

A  weakness  of  problem  (PU)  is  that  it  assumes  unlimited  capacity 
available  at  each  facility  in  terms  of  the  activities  using  a  given 
facility.  In  practice,  a  facility  may  not  have  the  capability  to 
serve*  every  activity,  and  may  have  restrictions  as  to  the  total 
capacity  available  to  handle  more  than  one  activity. 
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Lor  s  ”  capacity  available  at  farilitv  k  ,  atul 

d..^  -  capacity  required  at  facility  k  for  activity  j 
when  activity  j  uses  design  i  . 

If  design  i  does  not  include  facility  k  ,  then  d.  «  0  for  all  j 

k 

Define  the  decision  variable  y.  as: 

J  k 

y^  =  1  if  facility  k  is  used, 

=  0  otherwise  . 

Then  the  assignment  problem  [Gross,  Pinkus,  and  Soland  (1979)], 

called  problem  (P)  is  to  find  x..  and  y.  values  that: 

S  ij  Jk 


Minimize  I  E  a#.  x..  +  E  b,  y, 

1=1  j=1  u  k=1 


subject  to 


E  x.  .  =  1 


j*l» • • . ,n 


1  z  x,^l  k=l,...,p 


i-1  J-l 


*ijk  ij  -  °k'k 


x..,  y.  =  0  or  1  for  all  i>j,k 
ij  k  >J 


Constraints  (5)  of  problem  (P)  ensure  that  the  capacities  available  at 
the  facilities  are  not  violated.  Problem  (P)  is,  thus,  a  raul t iactivity 
multifacility  capacity-constrained  0-1  assignment  problem,  as  compared 
to  problem  (PU)  which  is  uncapacitated.  In  problem  (P) ,  constraints  (2) 
along  with  the  part  of  constraints  (6)  involving  the  fs  ensure  that 

each  activity  is  assigned  to  a  single  design.  Of  course,  the  optimal 
solution  may  result  in  the  use  of  more  than  one  design. 

For  an  example  of  five  facilities  and  three  designs  as  shown 

in  Figures  1(b)  to  1(d),  and  four  activities;  the  matrix  [a.. lb.  Id...] 

ij  1  k1  ljk 

is  as  shown  in  Figure  3, 

1.3.1  Comparison  with  the  uncapacitated  assignment  problem. 

Comparison  of  the  capacitated  problem  (P)  with  the  uncapacitated 
problem  (PU)  shows  that  the  objective  functions  (1)  and  (4)  are  equiva¬ 
lent  and  constraints  (2)  in  each  are  the  same.  Constraints  (5)  serve 
to  impose  the  capacity  constraints  and  at  the  same  time,  for  a  given 
design,  the  relevant  facilities  are  forced  in  the  solution.  For  an 


t*.  i  >^v  • 
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x  equal  to  1  ,  ail  the  facilities  with  d .  . ,  >0  must  have  y, 

ij  ijk  Jk 

values  equal  to  1  in  order  to  satisfy  (5)  and  the  corresponding  fixed 

costs  t>k  are  therefore  included  in  (4).  If  y^  -  0  and  d  ^  >  ® 


then  x  .  . 

ij 


must  be  0  in  order  to  satisfy  (5). 


Problem  (P)  has  been  formulated  as  a  0-1  linear  programming 
problem  whereas  problem  (PU)  was  formulated  as  a  0-1  nonlinear 
programming  problem. 


Note  that  problem  (PU)  can  be  easily  obtained  as  a  special 

case  of  problem  (P)  by  letting  d..,  equal  1  (for  all  j  )  if 

ij  k. 

design  i  uses  facility  k  ,  and  setting  all  equal  to  n  .  In 


other  words,  the  corresponding  formulation  is  to  find 
values  that: 


x .  . 
ij 


and 


y 


k 


Minimiz  e 


m  n  p 

L  £  a .  .  x .  .  +  I  b.  y, 
ij  ij  ,  ,  k  k 


i=l  j=l 


k=l 


(PI) 


m 

t 

i=l 


subject  to  £  x  ,  =  1 
ij 


j  1 »  •  •  •  t  n 


m 


l  x  .  <  ny,  k=l, - p 


i=i  'ik  j-i  ' ij  “  k 


x..,y,  =  0  or  1  for  all  i,j,k 

ij*  k 


^where 


'ik 


1  if  design  i  uses  facility  k  , 
0  otherwise  . 


(4) 


(2) 

(7) 

(6) 


1,3.2  Comparison  with  the  fixed-charge  locat ion-allocat ion  problem. 

Problem  (P)  bears  a  resemblance  to  the  well-known  fixed-charge 
location-allocation  problem  or  capacitated  facility  location  problem 
[Geoff rion  (1975);  Ross  and  Soland  (1977)].  There  are,  however, 
very  significant  differences  between  the  two.  In  order  to  point  out 
these  differences,  here  is  a  statement  of  the  location-allocation 
problem  (IJ\)  as  given  by  Gross,  Pinkus,  and  Soland  (1979)  in  a  form 


_  10 


similar  to  that  of  problem  (P)  . 

Find  and  values  that 


(LA)  < 


Minimize 


subject  to 


P  n  p 

z  z  xki  +  z  bkyk 

k=i  j=i  J  J  k=i 


m 

Z  x 

k=l 


kj 


j=l, . 


j=l 


djxkji  Vk 


xkj  >  0  ,  yk  =  o  or  1  for  all  j 

.ind  k 


(8) 


(9) 


(10) 


(ID 


Here  x^.  represents  the  fraction  of  customer  (activity)  j's  demand 

that  is  supplied  by  a  facility  at  location  k. 

The  most  important  distinction  between  problem  (LA)  and  problem 
(P)  is  the  relationship  between  assignments  and  facilities.  In  problem 
(LA)  there  is  a  direct  connection  between  the  assignments  made  and  the 
facilities  required,  and  each  assignment  affects  cnly  one  facility,  i.e., 
the  assignment  x^_.  >  0  has  a  bearing  on  only  one  of  the  constraints 

(10)  and,  therefore,  on  only  one  variable  y^  .  On  the  other  hand, 

in  problem  (P),  the  connection  between  the  assignments  made  and  the 
facilities  required  is  indirect,  and  each  assignment  can  affect  several 
facilities,  i  e.,  the  assignment  x__  =  1  bears  on  all  of  the  constraints 

(5)  for  which  d.  >  0  and,  therefore,  on  several  variables  y,  . 

1 J  K.  K. 

Another  distinction  is  the  relative  difficulty  of  the  two  problems. 
While  problem  (LA)  is  not  easy  to  solve,  branch-and-bound  approaches 
have  been  successful  in  dealing  with  it  because  oncz  values  are  specified 
for  the  ,  the  x.^  are  found  by  solving  a  transportation  problem. 

Problem  (LA)  becomes  more  difficult  if  the  constraints  x,  .  >  0  in  (10) 

kj  - 

are  replaced  by  x  =  0  or  1  in  order  to  preclude  supply  of  customer 
kj 

(activity)  j?s  demand  by  more  than  one  facility.  With  this  change. 
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problem  (LA)  may  be  treated  as  a  generalized  assignment  problem  and 
is  soLvable  using  an  efficient  branch-and-bound  algorithm  [Ross  and 
Soland  (1977)].  Problem  (P)  is  more  difficult  than  this  variation 
of  problem  (LA)  because  of  the  above  stated  indirect  connection 
between  the  assignments  and  the  facilities.  Even  after  values  have 
been  specified  for  all  the  y^  ,  problem  (P)  remains  a  difficult 

0-1  linear  programming  problem  because  of  the  interaction  of  the 
const ra int s . 

1,3.3  Solving  problem  (P). 

The  capacitated  problem  (P)  has  mn+p  0-1  variables  and 

n+p  constraints,  so  the  problem  dimensions  may  be  large  from  practical 
considerat ions.  For  example,  with  m^n^O  and  p=20  ,  problem  (P) 
has  920  variables  and  50  constraints.  The  0-1  LP  computer  codes 
generally  available  are  limited  in  terms  of  problem  size.  For  example, 
the  code  used  by  Gross,  Pinkus,  and  Soland  (1979)  can  handle  up  to  40 
variables  and  20  constraints.  A  better  and  more  efficient  code 
[Geoff rion  and  Nelson  (1968)]  allows  up  to  90  variables  and  50 
constraints.  This  fact,  together  with  the  structure  of  problem  (P) 
suggests  that  a  specialized  algorithm  could  be  developed  that  would 
be  more  efficient  for  practical  problems  than  the  general  integer 
linear  programming  algorithms  (on  which  the  available  codes  are 
based) . 

With  the  above  background  in  mind,  the  development  of  the 
solution  algorltlim  and  the  computer  program  to  solve  problem  (P) 
was  undertaken  and  is  described  in  Chapters  2  through  4. 

1.4  Areas  of  Application 

The  solution  algorithm  and  the  computer  program  are  designed 
to  solve  a  mu  1 1 iuc t ivi ty  multifacility  capacity-constrained  0-1 
assignment  problem,  i.e.,  one  which  can  be  formulated  as  problem  (P) . 

The  basic  elements  of  such  a  problem  are  activities  that  must 
be  assigned,  facilities  and  their  meaningful  configurations 


i 


i 

t 


V 

'i 

! 


i 


«a 


represented  as  designs,  the  fixed  and  variable  costs,  and  the 
capacity  requirements  of  the  activities. 


The  formulation  (P)  applies  to  existing  and/or  proposed 
facilities.  In  other  words,  it  is  useful  for  a  situation  where  the 
decision  may  be  to  delete  some  of  the  existing  facilities,  as  well  as 
for  a  situation  where  the  decision  involves  a  selection  out  of  a  set 
of  proposed  facilities. 

Table  1  includes  examples  of  areas  where  formulation  (P)  is 
applicable.  Within  each  application  area,  activities  and  facilities 
are  identified.  The  implications  of  designs,  variable  costs,  and 
fixed  costs  are  apparent. 


Obtaining  the  values  of  the  data  elements  b,  ,  d  .  ,  ,  s  , 

k  ijk  k 

and  in  particular  a__  ,  can  be  a  simple  or  a  complex  exercise  depending 

on  the  particulars  of  the  application,  and  the  nature  of  the  components 
comprising  these  elements.  For  example,  in  designing  multi-echelon 
inventory  sy items  [Gross,  Pinkus,  and  Soland  (1979)],  a_  represents 

the  inventory  cost  of  product  (activity)  j  using  echelon  structure 
(design)  i  and  b^  represents  the  fixed  cost  of  installation  (facility) 

k  .  The  inventory  cost  a_  includes  the  cost  of  procurement,  carrying 

inventory,  filling  orders,  and  stockouts.  The  value  a^.  ,  and  associated 

inventory  stockage  policies,  are  arrived  at  by  solving  a  multi-echelon 
inventory  problem.  In  other  words,  for  product  j  stocked  under  echelon 


structure  i 
installat Ion 


,  cptimal  inventory  policies  are  determined,  at  each 


of  the  structure,  which  yields  a 


ij 


The  facility  fixed 


cost  b^  includes  the  capital  expenditure  for  building  the  installation, 


along  with  a  number  of  fixed  costs  associated  with  operating  it,  such  as 
administrative  expenses,  the  expense  of  renting  the  facility  (if  it  is 
not  built),  .nd  certain  other  fixed  operating  expenses. 


In  the  case  of  designing  a  support  system  for  repairable 
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items  [Gross  and  Pinkus  (1979)],  represents  the  total  variable 

cost  if  unit  type  (activity)  j  is  repaired  under  design  i  .  The 
set  of  parameters  taken  into  consideration  to  compute  this  cost  for 
each  unit  type  includes  such  things  as  varying  population  sizes, 
failure  rates,  average  repair  times,  costs  associated  with  their 
repair,  the  purchase  and  storage  of  spares,  the  purchase  of  repair 
channels,  and  travel  to  depots  (facilities)  for  repair.  A  computer 
program  is  used  to  solve  a  spares  and  server  provisioning  problem,  and 
the  results  provide  the  basic  information  to  compute  a_  . 


Thus,  in  general,  the  data  elements  of  problem  (P)  may  be  obtained 
directly  and/or  by  solving  other  related  problem(s);  it  depends  on  the 
definition  and  the  nature  of  the  components  comprising  these  data  elements 
for  a  specific  application  area. 


2 .  DEVELOPMENT  OF  THE  SOLUTION  ALGORITHM 


The  solution  algorithm  that  has  been  developed  tt  solve  problem  (P) 
is  a  braneh-and-bound  procedure  which  makes  use  of  Lagrangian  relaxation 
as  a  basic  step. 

This  chapter  considers  two  different  Lagrangian  relaxations  of 
problem  (P),  their  general  characteristics,  and  some  useful  results  leading 
to  the  specific  case  of  Lagrangian  relaxation  util: zed  in  the  solution 
algor  i  tlim. 

2 . 1  Lagrangian  Relaxation 

Taking  a  set  of  '’complicating'’  constraints  of  a  general  mixed- 
integer  program  into  the  objective  function  in  a  Lagrangian  fashion  (with 
fixed  multipliers)  results  in  a  "Lagrangian  relaxation"  of  the  original 
problem  [Geoff rion  (1974)].  The  relaxed  problem  is  easy  to  solve  compared 
to  the  original  problem,  and  provides  a  lower  bound  (for  minimization 
problems)  on  the  optimal  value  of  the  original  problem. 

Although  the  use  of  Lagrangian  relaxation  in  discrete  optimization 
has  been  reported  prior  to  1970  [e.g.,  Lorie  and  Savage  { 1955),  Everett 
(1963),  and  Gilmore  and  Gomory  (1963)],  the  "birth"  of  the  Lagrangian 
approach  as  it  exists  today  [Fisher  (1978)]  occurred  in  1970  with  the 
successful  application  of  Lagrangian  relaxations  to  the  traveling  salesman 
problem  [Held  and  Karp  (1970,  1971)].  This  was  followed  by  application 
of  Lagrangian  relaxation  to  scheduling  problems  [Fisher  and  Schrage  (1972),  and 
Fisher  (1973,  1976)],  the  general  integer  programming  problem  [Shapiro 
(1971),  and  Fisher  and  Shapiro  (1974)]  and  the  generalized  assignment  problem 
[Ross  and  Soland  (1975)].  Table  2  lists  the  applications  of  Lagrangian 
relaxation  as  given  by  Fisher  (1978).  A  review  of  Lagrangian  relaxation 
is  also  provided  by  Shapiro  (1977)  and  Christofides  (1980). 
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TABLE  2 

* 

APPLICATIONS  OF  LAGRANGIAN  RELAXATION 


Problem 

Researchers 

Lagrangjan  Problem 

TRAVELING  SALESMAN 

Symme tr ic 

Held  &  Karp  (1970,  1971) 

Spanning  Tree 

Asymmetric 

Bazarra  &  Goode  (1977) 

Spanning  Tree 

Symmetric 

Balas  &  Christofides  (1976) 

Perfect  2-Matching 

Asymmetric 

Balas  &  Christofides  (1976) 

Assignment 

SCHEDULING 

n|m  Weighted 

Tardiness 

Fisher  (1973) 

Pseudo-Polynomial 

1  Machine  Weight 
Tardiness 

Fisher  (1976) 

Dynamic  Programming 

Pseudo-Polynomial  DP 

Power  Generation 

Systems 

Muckstadt  &  Koenig  (1977) 

Pseudo-Polynomial  DP 

GENERAL  IP 

Unbounded  Variables 

Fisher  &  Shapiro  (1974) 

Group  Problem 

Unbounded  Variables 

Burdet  &  Johnson  (1976) 

Group  Problem 

0-1  Variables 

Etcheberry,  et.  al.  (1978) 

0  -  1  GUB 

LOCATION 

Uncapacitated 

Cornuejols,  Fisher,  & 

Nemhauser  (1977) 

0-1  VUB 

Capacitated 

Geoffrion  &  McBride  (1977) 

0  -  1  VUB 

Databases  in 

Computer  Networks 

Fisher  &  Kochbaum  (1978) 

0  -  1  VUB 

GENERALIZED  ASSIGNMENT 

Ross  &  Soland  (1975) 

Knapsack 

Chalmet  &  Gelders  (1976) 

Knapsack,  0-1  GUB 

SET  COVERING-PARTITIONING 

Cover ing 

Etcheberry  (1977) 

0  -  1  GUB 

Partitioning 

Nemhauser  &  Weber  (1978) 

Matching 

*Source:  Fisher  (1978) 
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4.1.1  Relaxing  Problem  (P) 

liv  dividing  constraints  (5)  by  s  and  letting  r  *  *  =  d  -  •  *,  /  c  > 

k  l  J  1 J  K '  s  r 

problem  (P)  can  be  restated  as  follows. 


Minimize  Z  Z  a.,  x..  +  Z  b,  y, 

i-l  j=l  k=1  k  k 


sub j  ec  t  to  Z  x . .  =  1 
i-i 


j ,n 


£  ^  rm  \s  -  yk  k*1 . p 


x.,»  y .  =  0  or  1  for  all  i,j,k 

ij  k 


A  Lagrangian  relaxation  (LR^)  of  problem  (P)  relative  to  constraints 


(2)  is  obtained  as 


Minimize  Z  Z  a . .  x . .  +  Z 
1=1  j  ~  1  1J  [<;  = 


(LRu>l  subject  to  ^^rijk  x..  <  yR 


b  y.  -  E  u.  (  I  x.  .  -  l\ 

1  k  k  j-1  3  \i-i  13  / 


x .  .  ,  y,  -  0  or  1 
lj  k 


k=l, . . . ,p 


for  all  i , j , k 


where  the  u  are  Lagrange  multipliers;  it  follows  that  the 

optimal  value  of  problem  (LR^)  is  a  Zower  bound  on  the  optimal  value  of 

problem  (P) ,  i.e.,  Z(LR  )  <  Z(P)  .  We  will  continue  to  use  this  notation 

u  — 


in  which  /.(•)  is  the  optimal  value  of  problem  (•) 


Another  Lagrangian  relaxation  (LR^)  of  problem  (P) ,  relative  to 


constraints  (5’),  is  obtained  as 
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2.1.2  General  Charac teristics 

A  Lagrangian  relaxation  provides  a  lower  bound  on  the  optimal  value 
of  the  orginal  problem,  i.e.,  in  our  ease  Z(LR^)  Z(P)  and  Z(LR^)  <_  Z(P) 
The  usefulness  of  a  Lagrangian  relaxation  depends  du  the  closeness  of  this 
lower  bound  to  the  optimal  value  of  the  orginal  problem.  However,  the 
relaxation  must  be  "easy”  to  solve  relative  to  the  original  problem.  We 
observe  that  the  optimal  value  of  in  problem  (LK^  is  1  if 

(v,  -  h  )  >  0  and  0  if  (v,  -  b,  )  <  0  ,  and  then  prol lem  (LR  )  reduces 
k  k  —  k  k  —  v 

to  n  0-1  "multiple  choice"  problems  which  are  very  easy  to  solve. 

On  the  other  hand,  problem  (LR^)  reduces  to  k  0-1  knapsack  problems. 

However,  these  problems  are  not  independent  because  of  the  interaction 

of  constraints  (51)  and  the  indirect  relationship  described  earlier  in 

Section  1.3  between  the  assignments  and  the  facilities.  In  view  of  this 

complexity,  relaxation  (LR^)  will  not  be  considered  further. 

The  choice  of  Lagrange  multipliers  in  relaxation  (LR^)  should 

be  such  that  Z(LR^)  is  as  large  as  possible  and  hence  as  close  as 

possible  to  Z(P)  in  view  of  the  relationship  Z(Lk^)  _<  Z(P)  .  In  other 

words,  an  equivalent  problem  is  to  find  a  vector  v  f represent ing  v^,  v^* 


to 
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(») 


Maximize 


v  >  0 


Obviously,  Z(LR  )  <  Z(I))  <  Z(P)  . 


(K) 


The  general  properties  of  Lagrangian  relaxation  have  been  well 
described  in  the  literature  [e.g.,  Geoffrion  (1974),  Geoffrion  and  McBride 
(197b),  and  Fisher  (1980)].  Some  oi  these  properties  relating  the 
Lagrangian  relaxation  and  the  usual  LP  relaxation  are  stated  below. 

The  LP  relaxation  (P)  of  problem  (P)  is  obtained  by  relaxing  the 
integrality  constraints  (6),  i.e.,  the  formulation  (P)  is 

^  m  n  p 


Minimize  Z  Z  a  x  +  Z  b,  y.  (4) 

ij  ij  k  k 


i=l 

j-i  1J 

IJ 

k=l 

k  k 

m 

subject  to 

Z  x  .  * 
i=i  iJ 

1 

c 

H 

II 

(2) 

(!>)  * 

m 

n 

Z 

i=l 

E  r. 
j-i  *3k 

X  .  . 

IJ 

1  yk 

k=l,  .  . .  ,p 

(5’) 

k=l, .  . .  ,p 

(15) 

x .  .  , 
ij 

o 

A  | 

for  all  i, j  9 k 

(16) 

Note  that 

the  constraints 

X  .  . 

IJ 

<  1 

are  implicit  in 

constraints  (2) 

Also  consider  the  following  partial  convex  hull  relaxation  (P*) 
of  problem  (P) . 


Mini mize 

m 

Z 

i=l 

n  p 

Z  a  .  .  x  .  .  +  Z  b.  y 
i=1  iJ  iJ  k=1  k  k 

(*> 

([’*)  < 

subject  to 

m 

Z 

i=l 

n 

E  r  ,  . .  x . .  <  y.  k=l , . . . ,p 

j=1  ij  -  k 

(5’) 

k 

x  . 

1 

. ,  y.  £  convex  hull  {(2),  (6)} 

I  k 

(17) 
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Then  the  relationships  between  the  optimal  values  of  various 
problems  [Geoff rion  and  McBride  (1973) ]  are  as  follows. 

2(F)  Z(LR^)  <  max  Z(LR  )  =  Z(D)  «  Z(!>*)  <  Z(P)  (18) 

v  >  0 

'v  'v  %  % 

where  v  are  the  values  v^,  v2»*'*>Vp  a  dual  optimal  solution 
of  (P)  corresponding  to  constraints  (5*). 

Thus,  the  optimal  dual  solution  associated  with  the  usual  LP  relax¬ 
ation.  furnishes  a  choice  of  Lagrange  multipliers  such  that  the  associated 
Lagrangian  relaxation  is  at  least  as  tight  as  the  usual  LP  relaxation, 
and  generally  a  good  deal  tighter  and  even  as  tight  as  the  partial  convex 
hull  relaxation. 

Since  Z(D)  =  Z(P*)  ,  the  quality  of  the  bound  obtained  from  the 

Lagrangian  relaxation  depends  on  where  Z(P*)  lies  in  the  range  between 

Z(P)  and  Z(P)  .  It  turns  out  that  problem  (LR  )  possesses  the 

v 

"integrality  property, 11  i.e.,  the  optimal  value  of  problem  (LR^)  is  not 

altered  by  dropping  the  integrality  conditions  on  its  variables  and 
therefore  [Geoffrion  (1974)] 

Z(D)  =  Z(P*)  =  Z(P)  (19) 

Thus,  the  Lagrangian  relaxation  (LR^)  is  no  better  than  the  LP 

relaxation  (P) .  On  the  other  hand,  Lagrangian  relaxation  (LR  )  does 

not  possess  the  integrality  property  and,  hence,  could  provide  an  equal 
or  better  bound  than  the  LP  relaxation  (P)  ;  but  the  compu tational 
difficulties  do  not  favor  pursuing  formulation  (LR^)  . 

It  is  possible  to  consider  alternative  formulacions  of  problem 
(P)  with,  the  objective  of  obtaining  tighter  bounds.  This  aspect  is 
discussed  in  Chapter  6. 
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2.2  Sumo  Result.: 


We  now  turn  to  the  basic  question  of  choosing  Lagrange  multipliers 

v  so  that  ( LR  )  Is  optimal  to  the  extent  possible,  which  is  equivalent 

to  solving  problem  (D)  .  We  also  need  to  consider  this  question  when  5;ome 

of  the  x..  and  y,  variables  have  been  assigned  values  of  1  or  0,  i.e., 
1 J  R 

at  a  node  other  than  the  starting  or  '’root1'  node  in  the  fcranch-and-bound 
tree.  For  this  purpose,  some  terminology  is  defined  and  formulations 
corresponding  to  problems  (P) ,  (LK  )  and  (D)  are  first  developed.  Then 

some  important  results  pertaining  to  the  choice  of  Lagrange  multipliers 
will  be  proved.  Gavish  (1978)  provides  a  method  of  obtaining  the  ’best’ 
multipliers,  based  on  solving  an  equivalent  linear  programming  problem. 
Such  a  formulation  is  difficult  in  our  case,  and,  besides,  we  propose 
to  avoid  solving  LP  problems  in  our  branch-and -bound  procedure. 

Define  the  sets 


S  =  {(i»j)jx^j  has  an  assigned  value  of  1  or  0}  ,  and 

T  =  (k|y  has  an  assigned  value  of  1  or  0)  . 

These  sets  represent  the  partial  solution  of  problem  (?)  and  the  variables 
contained  in  these  sets  are  termed  f ixed  variables.  [Geoff rion  (1967)]. 
Let  S  and  T  represent  the  corresponding  complementary  sets,  i.e., 
comprised  o f  the  x.j  and  y^  variables,  which  have  not  been  assigned 

spec  11 ic  values  ; .ml,  t heref ore,  are  called  free  variables.  A  completion 

of  a  partial  solution  is  defined  as  a  solution  that  is  determined  by 

S  and  T  together  with  a  binary  specification  (0  or  1)  of  the  values 

oi  the  free  x..  and  y.  variables  from  sets  S  and  T  . 
ij  ^k 

Let  SUS  =  S  and  TVJT  *  Tj  . 

’onsider  a  partial  solution  to  problem  (P)  in  which  specific  values 
(of  1  or  0)  are  assigned  to  some  of  the 


x..  and  y.  such  that 
ij  y  k 
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m 

Z  x,  .  <_  1  Vj  > 

i=l  J 
(i,j)eS 


and 


Z  Z  r 

i  j 

(i> j) £S 


ijk 


x. . 
iJ 


-  'k 


VkeT 


Z  Z  r . ..  x. .  <  1  VkeT 

i  j  ijk  ij  - 
(i,  j)eS 


ana  such  that  x.  .  -  1  and  e  —  1  imply  tnat  keT  and  y,  “  1* 

ij  ik  k 

Recall  that,  by  definition,  eik  -  1  if  design  i  uses  facility  k  , 

and  e  ,  -  0  otherwise, 

ik 


The  problem  of  finding  an  optimal  completion  of  the  partial 
solution  of  problem  (P)  can  be  stated  as  follows. 


Minimize 


subject  to 


i  ]  a«  + 1  Vk +  55  + 1  Vk  (20) 

(i,j)eS 


k 

kef 


i  j 
(i,j)es 


k 

keT 


<h>  * 


E  x .  .  =  1  -  £  x 
i  i  U 

(i, j)eS  (i, j)eS 


Vj 


l  ]  rijk  xij  1  yk  '  l  J  rijk  xij  Vk 


U,J)eS 


(i.J)eS 


x,  .  ,  y  =  0  or  1  V(i,j)£S  ,  keT 

ij  k 


(21) 


(22) 


(23) 


We  call  this  problem  (P^)  where  H  indicates  the  node  in  the  branch-and- 
bound  tree. 


A  Lagrangian  relaxation  of  problem  (P^)  with  respect  to  constraints 
(22)  is  obtained  by  introducing  non-negative  Lagrange  multipliers  , 

k»l,2,...,p  ;  the  relaxation  is  then 
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Minimize  X  Y.  a  .  x  +  £  by,  +  E  li  a  x  +  Z  b  v 

i  j  _  u  ij  k  k  k  .  .  ij  ij  k  kyk 


(i, j)eS 


keT 


T-  v  y  -  n  r  . . ,  x  . .  -  Z  Z  r  x . . 

i  k  k  .  .  ljk  ij  .  .  l  ]  k  ij 

k  1  j  i  j 


(i. j)£S 


(i, j)£S 


suoiect  to  Z  x..  =  1  -  Z  x.. 

ij  .  ij 

1  i  J 

(i,j)£S  (i,j)ES 


x  .  . ,  y  =  0  or  1 
ij  k 


V(i,i)f.S  ,  keT 


Rearranging  (24),  and  using  the  relationship  =  TUT,  we  have  problem 


(LR.  ) 
v 


I  Minimize 


(LR  J 


l  Z  X . .  /a . .  +  I  V,  r .  .  \  +  Z  Z  x . .  ( a . .  +  Z  v,  r  . . ,  | 
i  j  _lJ  I  1J  kcT,k  «kl  i  j  Vl  ij  kET  k  Uk) 

[i,j )es  '  1  /  (ij)es  1 


k  (\  -  \)  ^  (\  -  \) 


subject  to  Z  x..  =  1  -  Z  x.. 

t  U  i  U 


(i, j)cS 


(1 » j )es 


x .  . ,  y  =  0  or  1 
ij  Jk 


V(i,j)eS  ,  keT 


Then  we  have  Z(LR^  _<  Z(Pj)  •  An  important  problem  is  the  choice  of 
Lagrange  multipliers  v^,v^,...,  ,  represented  by  vector  v  ,  that 

maximize  Z(LK^  y)  ,  i.e.,  the  problem  (D^)  : 


Maxi  ini?; 


“  [Z"-RCv>] 


v  >  0 
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We  now  state  and  prove  some  theorems  related  to  the  choice  of  Lagrange 
multipliers  V|,V2’’**,Vp  * 

Theorem  1;  There  exists  an  optimal  solution  to  problem  (D)  in 
which  _>  for  all  k  . 

Proof :  Suppose  ,  in  an  optimal  solution  to  problem 

(D),  i.e.,  Z (D)  =  Z(LR  )  where  /  *  <  b  . 


Recall  that 

Z  (LR  .  )  *  Min  E  I  x,,  (a.,  +  E  vf  r,,.  U  «  y.  /v*  -  b.  ^ 

v*  i  j  V 1J  k  k  13k)  k  k  [ k  v 

s.t.  Zx..=l  Yj  (2) 

i  iJ 

x,,,  y  -  0  or  1  Vi,  j  * k  (6) 

lj  K 

For  v*  <  b^  ,  the  optimal  value  of  y^  is  0  ,  and  the  term  -  (v*  -  b^) 

in  the  objective  function  is  0. 

Consider  what  happens  if  we  increase  v*  to  .  Call  the 

resulting  vector  .  Consider  problem  (LR^)  .  The  optimal  value  of 

y^  in  problem  (LR^)  is  0  or  1  ,  and  the  term  -  y^  (v^  -  b^)  is  0  . 

However,  the  optimal  value  of  y^  is  the  same  in  problems  (LR^*)  and 

(LR  )  for  all  k  >  1  .  Therefore,  the  quantity  L  y^  (v  -  b  )  is  the 
—  k 

same  at  the  optimal  solution  for  both  v  =  v*  and  v  =  v  . 

Since  >  v*  ,  we  note  that  in  the  objective  function, 

P  P 

a  +  2;  V  r  >  a.,  +  £  vf  r 
lj  k=1  ~ *  ijk  -  lj  k=1  k  ijk 


and  therefore  Z(LRy)  >  Z(LRy*)  . 


Vi,j  , 
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It  follows  that  there  is  an  optimal  solution  to  problem  (D) 
in  which  . 

Since  the  choice  of  k=l  was  arbitrary,  the  same  result  holds  for 

any  value  of  k,  k=l,*..,p  ;  hence,  there  exists  an  optimal  solution  to 

problem  (D)  in  which  v,  >b  for  all  k. 

k  —  k 

Theorem _ 2:  There  exists  an  optimal  solution  to  problem  (D,0  in 

which  _>  b^  if  (i)  keT  or  (ii)  keT  and  =  0  . 


Proof :  Suppose  in  an  optimal  solution  to  problem 

(D  )  ,  i.e.,  Z(D. )  =  Z(LR„  *)  where  v*  <  b.  . 
x,  i  £,v  1  1 

Then  k=l  can  be  such  that  ktT  or  keT  . 

Case  (i) :  Let  ksT  . 


term 


2  6 


Lot  v*  bo  increased  to  ;  call  the  resulting  vector  v  . 

Consider  problem  (LR  )  .  The  optimal  value  of  y„  in  (LRn  )  is 

£ ,  v  1  £,v 

0  or  1  ,  then  the  term  -  y  (v^  -  b^  is  0  .  For  k  >  1  ,  the 

optimal  value  of  y  being  the  same  in  (LR.  *)  and  (LR.  )  ,  we 
k  £,v  £,v_ 

find  that  Z  y  (v  -  b  )  is  the  same  at  the  optimal  solution  for  both 
keT  K  k  k 


and  v  = 


But  v_-^  >  v^  ;  therefore 


a  .  .  +  Z  v.  r  .  >  a  .  .  +  Z  v,*  r .  .  _ 

1J  ljk~  1J  keTl  k  ljk 


V(i,j)eS  and  (i,j)eS 


Hence,  Z(LR  )  Z(LR„  *)  ,  wherefrom  it  follows  that  there  is  an 
optimal  solution  to  (D,)  in  which  v  >  b  .  Sr.nce  the  choice  of 

A'  1  "  1 

k=l  was  arbitrary,  the  same  results  hold  for  any  vaJue  of  k,  kef  . 
Hence,  there  exists  an  optimal  solution  to  problem  (0^)  in  which 

v,  >  b,  for  ail  kef  . 
k.  —  k 

Case  (ii):  Let  keT  and  y  =  0 


Considering  problem  (LR0  *)  ,  for  k  =  1  ,  <  b.  ,  and 

v  11 

y^  =  0  ,  the  term  -  y^(v*  “  in  the  objective  function  is  0  . 

Increase  v^  to  and  call  the  resulting  vector  v_  .  The 

term  -  y^(v ^  -  b^)  is  0  .  For  k  >  1  ,  the  optimal  values  of  y^ 


are  the  same  in  problems  (LR.  *)  and  (LRP  )  .  Therefore  Z 

*’V  kcT 


y^(v^  -  b.^)  is  the  same  at  the  optimal  solution  for  both 


v  =  v 


and  v  -  v 


Since  v.,  >  v 


1  » 


aii  +  Y-  ^kriik  -  Uii  +  E  Vk  riik  VOJKS  and  (i, j)eS 
J  ktT1  k  Jk  J  kcTL  k  Jk 
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Therefore  Z(LR,  )  >  Z(LR.  *)  .  It  follows  that  there  exists  an  optimal 

•*-  ,  V  —  X,  ,  V 

solution  to  (D^)  in  which  v^  >  .  The  choice  of  k=l  being  arbitrary, 

the  same  results  hold  for  any  value  of  k  ,  kcT  and  -  0  ;  which  proves 

case  (ii)  of  the  Theorem. 

It  may  be  added  that  there  is  another  possibility  which  complements 
case  (ii)  of  Theorem  2,  i.e.,  if  kcT  and  y  =  I  .  We  treat  this  possi¬ 
bility  as  a  conjecture  since  a  result  similar  to  the  one  above  could  not 
be  proved,  as  discussed  now. 

With  y^  =  1  and  v*  b^  ,  uTe  observe  from  problem  (LR,„  v*) 
that  for  a  solution  vector  X  '  (with  elements  x_*)  anc^  Y*  (with 
elements  y  v* | Yi 1  and  y0*,...,y  *  s  0  or  1)  , 

J 1  -p‘^1  l  J  p 


i  J 

(i. j)eS 


Z(LR,  *)  =  Z  E  x.  .*  [  a.  .  +  v*  r  .  ..  +  Z  v  *  r 


;  ij  l  ij  1  ijl 


k  ij  k 


i  J 

(i,j)eS 


4*  Z  E  x  .  .  *  [  a  .  .  +  v  *  r  .  . ,  +  E  v,  *  r  .  . 


ij  V  ij  1  ij  i 


k  ijk 


-  £  y 

kCT 


(*£  -  \)  ~  -  \)  V  ( 


'  v  *  -  b 
k  k  k 


Since  v^*  <  b^  and  y^  =  i  the  term  -  y^(v^*  -  b^)  is  positive. 

If  we  raise  v^*  to  ,  say  ,  the  term  -  y^(v^  -  b^)  is  0  . 

The  difference  between  Z(LR  )  and  the  objective  function 

*  * 

value  of  problem  (LR^  with  X=X  and  Y=Y  is 

=»  E  E  x.  .  *  v*  r  .  +  (b  -  v  *)  -  X  E  x..*  b  r 

.  .  ij  J  i  j  1  1  1  .  .  ij  ij  1 

i  j  J  J  i  J 

=  (b  -  V  *)  -  Z  l  X  *  (bx  -  vx*)  r  . 

1  J 

This  difference  can  be  either  negative  or  positive,  and  s)  we  cannot 
conclude  that  there  is  an  optimal  solution  to  problem  (D^)  in  which 
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V1  >  ‘J1  *  kV  believe  Lhis  conclusion  to  be  false. 

it  it 

Theorem  3:  Let  (X  ,  Y  )  solve  problem  (LR^)  f°r 

/  *  * 

for  all  k  .  If  (X  ,  Y  )  is  feasible  for  problem 

(P) ,  there  exists  an  optimal  solution  to  problem  (D) 

in  which  v,  =  b.  for  all  k  . 
k  k 

Proof :  In  view  of  Theorem  1,  there  exisLs  an  optimal  solution 

to  (D)  in  which  v^  >_  for  all  k  *  i.e.,  v  >  b  . 

Let  v  be  such  an  optimal  v  .  We  will  show  that 
Z(LRv)  Z(LR^)  ,  from  which  it  follows  that  v  -  b 

solves  problem  (D)  . 

Recall  that 

Z(LR  )  =  Min  E  E  x  (a  +  E  v,  r,  i  -  I  v,  /v,  -  b.  \ 

Z  X,Y  i  j  ^  V  ij  k  “*  ljk  /  k  k  V_k  k) 


s .  t ,  Ex..  1 
,  ij 


ij 


yk 


=  0  or  1 


Vj 

Vi, j,k 


(2) 

(6) 


Since  v>b,  y  =  1  Vk  is  an  optimal  choice. 
-  k 


Hence,  Z(LR  )  -  Min  E  E  x., 

v  -  -■ 

-  X 


E  E  x  .  .  / a 

ij ij  \ 


+  E  v,  r 


iJ  k  “k  ijic 


s.  t.  Ex  =1 

iJ 


VJ 


x  =  0  or  i  Yi,j 

ij 


(2) 

(6a) 


Now  consider  (LR^)  •  Since  v  =  b  ,  the  last  term  of  the  objective 
function  drops  out,  and  we  have 


Z(LR,  )  =  Min 
b  X 


in  E  E  x  fa 
,Y  ij  ^  \i 


+  E  b,  r 
ij  k  k  ljk 


) 


subject  to  (2)  and  (6) 


Min 

X 
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Z  X  x  . 

(a  .  +  I  h  r  ,  J 

i  j  iJ 

\  U  k  k  Uk/ 

subj  ect 

to  (2)  and  (6a) 

x.  .*  fa 

+  Z  b  r  J 

1J  \ 

ij  k  k  ijk/ 

where  X  with  elements  x..*  is  the  minimizing  solution  vector 

ij 


which  satisfies  (2)  and  (6a)  . 


*  * 

Now  (X  ,  Y  )  feasible  for  (P)  implies  that 


l  Z  x .  . x  r  ,  t  <  v.  *  <  1  V  k 

xi  ij  k  —  '  k  — 

x  J  J  J 


Hence,  Z 
k 

(*  • bk) 

Z  Z  X. ,*  r 

.  .  ij  ijk 

l  j  J 

or  Z  | 

k  ' 

Z  Z  x .  .  *  r  .  .  * 

i  j  U  ijk 

Rewriting, 

Z  ( LR  )  = 

V 

Min  ZZx. 

X  ij  i] 

’-k  ~  bk 


a  .  +  Z  r  .  . ,  /b.  +  / 

lj  k  \k  \- 


subject  to  (2)  and  (6a) 


=  Min 
X 


in  I  I  x .  .  /a .  .  +  Z  r  .  f  b,  \ 

k  j  i  j  ^  k  i:)k  7 

+  \  Z  Z  x .  .  r  .  . .  -  Z  /  v 

ki  t  i « 'Jk  k  (-k 


-  b. 


subject  to  (2)  and  (6a) 


=  ZZx..*  fa .  .  4-  Z  r  b.  \ 
i  j  U  (  Xj  k  Ijk  k) 

+  1  A  -  b.  \t.  Z  X,  .*  r.  ..  -  I  /v 

k  rk  k/ij  ij  i-Jk  k  [~ 

<  Z  I  X.  .*  /  a . .  +  Z  r  .  . ,  b,  \  =  Z(LR,  ) 

~  i  j  1J  (  1J  k  ljk  k/  % 


'A  ****,♦.  *  *  .  #  -  Ate  "•?**** . 
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by  (27),  or  ZfLR^)  Z(LR^)  ;  it  follows  that  v  =  b  solves 
problem  (D) . 

2 . 3  Relaxation  (PR  ) 

Theorems  1  and  3  are  useful  in  providing  a  choice  of  Lagrange 

multipliers  as  a  starting  point  in  solving  a  relaxation  of  problem  (P) 

at  the  root  node.  Theorem  2,  similar  to  Theorem  1,  provides  results  for 

a  partial  solution  of  problem  (P) ,  i.e.,  at  a  node  other  than  the  root 

node  where  some  of  the  x..  and  y,  have  been  fixed  at  1  or  0  . 

ij  k 

Theorem  1  is  important  in  pointing  out  that  a  certain  set  of  Lagrange 
multipliers  v  such  that  v.  >  b,  for  all  k  would  provide  an  optimal 

K  K. 

choice.  Theorem  3  narrows  this  choice  to  v  =  b  for  all  k  for  a 

k  <. 

specific  situation,  i.e.,  when  the  resulting  solution  is  feasible  for 
problem  (?) . 

Letting  -  b^  for  all  k  ,  problem  (LRv)  becomes: 


r 


Minimize 

XL  c  .  . 

i  j  1J 

x .  . 
ij 

(28) 

r* 

50 

cr 

^ _ 

subject  to 

E  x.  .  »  1 
i  1J 

Vi 

(2) 

x. .  -  0 
ij 

or  1 

Vi.j 

(6a) 

where  c  .  .  = 
iJ 

a .  .  +  E 
k 

b.  r  .  ..  . 

x  ijk 

(29) 

Note  that  problem  (LR. 

b 

)  is  very 

easy  to  solve;  its  optimal  value 

is  just 

the  sum 

of  the  minimum 

(over  i) 

o 

p* 

t — i. 

o 

t-1 

all  j  ,  i.e., 

=  E  min 

j  i 

(c  .  .} 

(30) 

We  solve  this  problem  as  a  starting  point  at  the  root  node  in  our  branch- 
and-bound  procedure. 
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5ki  Max  {aR!l  ^  } 


so  y  is  fixed  at  i  and  kcT  if  6^  =  1 


To  return  to  the  relaxation,  problem  (LR.  ),  we  must  make  a  choice 

X/  ,  v 

of  the  vector  v  of  Lagrange  multipliers.  Of  course,  we  would  like  to 
use  an  optimal  choice,  i.e.,  a  vector  v  that  solves  problem  (D^)  • 

Recall,  however,  that  Theorem  2  did  not  provide  us  any  useful  information 
about  the  optimal  value  of  v^  if  keT  and  ”  1  .  To  simplify  our 

approacn  and  have  recourse  to  the  results  of  Theorems  1  and  3,  we  choose 
=  0  if  keT  and  y^  -  1  .  Note  that  there  a ~e  no  kcT  such  that 

V  =  0  because  of  practical  considerations  and  because  our  branching  rule 
Jk 

only  results  in  fixing  y^  values  at  1  .  Problem  (LR^  v)  now  takes  the 


Minimise 


Z  Z  x .  .  /  a  . 

m  1J  \ 


+KcfVk  ^  h  -  ^V”’ 


(LR^  -)  i  subject  to  Z  x_  =  1  , 


x.  .,  y,  a  0  or  1  for  all  (i,j)e?  ,  keT  . 

1  j  K. 


Note  that  in  this  problem  (LR0  -) ,  v  =  0  if  keT  .  Also  note  how  closely 

it  resembles  problem  (LR^),  the  relaxation  at  the  root  node.  As  in  that 

case,  we  would  like  the  lower  bound  ^(LR.  -)  to  be  as  large  as  possible, 

x. ,  v 


i. e. ,  we  seek  v  to 


Maximize 


v  >  () 


Because  of  the  close  similarity  of  problems  (LR^  -)  and  (LRv), 
it  is  possible  to  obtain  results  about  problem  (Dr)  that  are  analogous 


to  those  obtained  about  problem  (D) .  We  state  these  results  as  Theorems 
4  and  3.  Their  proofs  are  omitted  because  they  follow  precisely  the 
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This  specific  relaxation,  problem  (PR^) ,  is  ot  the  same  form  as 

problem  (LR,  )  and  is  equally  easy  to  solve  in  one  pass.  Its  optimal 

value  Z(PR^)  serves  as  the  lower  bound  at  node  S  .  Note  that  for  £=1  , 

problem  (PR0)  is  the  same  as  problem  (LR,  )  . 

x,  b 

It  is  clear  that  setting  each  Lagrange  multiplier  to 

for  kef  and  to  0  for  keT  is  not  generally  optimal  in  terms  of 
achieving  the  tightest  lower  bound  (except  as  per  Theorem  3).  But  it 
provides  a  good  starting  point  in  seeking  an  optimal  vector  v  and  it 
provides  an  easily  calculated  lower  bound  at  each  node  of  our  branch- 
and-bound  procedure.  The  question  of  how  to  improve  upon  this  choice 
of  multiplier  values  will  be  discussed  in  Chapter  6. 


3 .  METHODOLOGY  FRAMEWORK 


The  branch-and-bound  procedure/methodol ogy  developed  to  solve 
problem  (P)  uses  Lagrangian  relaxation  (PR^)  as  a  basic  step.  The 

branching  rule  dictates  which  x,.  variable  to  branch  on  at  each  node. 

ij 

in  addition,  there  are  certain  rules  (e.g.,  the  capacity  rule  and  the 
bounding  rule)  which  contribute,  significantly,  in  improving  the  overall 
efficiency  of  the  procedure. 


Some  basic  terms  such  as  fixed  and  free  variables,  partial  solution 
and  its  completion  were  introduced  in  the  previous  chapter.  This  chapter 
first  provides  a  preliminary  discussion  of  the  branch-and-bound  methodology, 
[Ceoffrion  (1967),  and  Geoff rion  and  Ilarsten  (1972)].  Representation  and 
storage  of  the  x..  variables  for  branching  and  backtracking  is  described 

in  order  to  pro\ ide  continuity  and  consistency  with  the  computer  program 
covered  in  Chapter  4.  This  is  followed  by  a  description  of  the  major 
components  of  the  branch-and-bound  methodology. 


Branching  and  backtracking  is  done  on  the  x..  variables.  The 

ij 

branching  commences  by  t ixing  the  variable  (selected  by  the  branching 

rule)  to  1  and  moving  to  the  left  branch  node.  When  backtracking,  we 
fix  the  corresponding  variable  at  0  and  move  to  the  right  branch 

nude  (if  the  right  branch  node  has  not  already  been  explored).  An 

variable  can  also  be  fixed  at  0  or  1  by  rules  other  than  the  branching 
rule.  The  capacity  rule  and  the  bounding  rule  are  two  such  rules  employed 
in  our  methodology. 


Figure  4a  shows  a  branch-and-hound  tree.  The  x  variables  fixed 

at  0  or  1  at  any  node  due  to  rules  other  than  the  branching  rule  are 
shown  in  parenthesis  at  the  appropriate  node. 
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,  Q  is  the  root  node  and  also  the  parent  node  for  nodes 
0  and  © 

>  ©  is  the  parent  node  for  nodes  ©  and  ©,  etc. 


Figure  4a. 

A  branch-and-bound  tree  illustration 


Node  (ft) 


Partial  Solution  (S 


1 

2 

3 

4 

5 

6 
7 


<P 

{103 v  -  301,  -  401} 

{103,  -  301,  -  401,  402} 

(103,  -  301,  -  401,  -  402  -  201,  101) 

(103,  -  301,  -  401,  -  402,  -  201,  101,  204) 

{103,  -  301,  -  401,  402,  -  201,  101.,-204) 


{-  103} 


Partial  solutions  for 


Figure  4b. 

the  above  illustration  (Figure  4 a) 
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For  problem  (PR^)  ,  a  partial  solution  corresponding  to  set  S  at 

node  l  ,  i.e  ,  S.,  contains  x..  variables  assigned  values  of  1  or 

£  ij 

0  .  For  simplicity  in  the  computer  program,  an  variable  fixed  at 

1  is  represented  as  (100  i  +  j)  ,  whereas  an  x  variable  fixed  at 

0  as  -  (100  1  +  j)  ,  e.g.,  x  9  *  1  and  x32  -  0  a:*e  represented  as 

302  and  -  302  respectively.  Since  branching  is  done  on  x_  variables, 

it  is  necessary  to  make  a  distinction  between  x_  variables  fixed  at  1 

due  to  the  branching  rule  and  those  fixed  at  1  due  to  the  other  rules. 

We  make  this  distinction  by  underlining  the  posit ive  number  to  represent 

an  x..  fixed  at  1  due  to  the  other  rul es .  For  example,  204,  -  301, 
ij 

_103  represent,  respectively,  x2^  =  1  due  to  the  branching  rule,  =  0 

due  to  the  branching  rule  or  any  other  rule,  and  -  1  due  to  a  rule 

other  than  the  branching  rule. 

Figure  4b  shows  the  partial  solutions  of  the  branch-and-bound 

tree  in  Figure  4a. 


Implicit  enumeration  involves  generating  a  sequence  of  partial 
solutions  and  simultaneously  considering  all  completions  of  each.  For 
our  minimization  problem,  we  start  with  an  initial  solution  having  a 
very  large  value  (Infinity)  as  an  initial  upper  bound.  As  the  computations 
proceed,  feasible  solutions  (those  satisfying  the  capacity  constraints) 
are  discovered  from  time  to  time,  and  the  best  one  yet  found  is  retained 
as  an  incumbent  solution  with  the  corresponding  value  as  the  best  upper 
bound .  It  may  happen  that  for  a  given  partial  solution  we  can 

determine  a  best  completion  of  ,  i.e.,  a  feasible  completion  that 

minimizes  the  objective  function  value  among  all  feasible  completions 
oi  .  If  such  a  best  feasible  completion  is  better  than  the  best 

upper  bound,  then  it  replaces  the  latter.  Or  we  may  be  able  to  determine 
that  lias  no  feasible  completion  better  than  the  incumbent.  In  either 

case,  we  can  f a t horn  .  (Various  situations  of  fathoming  and  back- 
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tracking  in  our  branch-and-bound  procedure  are  described  in  the 
following  discussion.)  All  completions  of  a  fathomed  partial  solution 
have  been  implicitly  enumerated  in  the  sense  that  they  can  be 

excluded  from  further  consideration  (with  the  exception  of  the  relevant 
best  feasible  solution  of  if  it  has  been  retained  as  the  best 

upper  bound) . 

In  our  branch-and-bound  procedure,  at  any  given  node  where  we  can 
fathom  ,  we  backtrack  to  the  parent  node  and  move  to  the  right-hand 

branch  (if  that  branch  has  not  already  been  explored)  by  fixing  the 
appropriate  x^  variable  at  0  .  However,  if  the  right-hand  branch 

has  already  been  explored,  we  continue  backtracking  to  a  parent  node  where 
we  can  move  to  a  right-hand  branch.  For  example,  in  Figure  4a,  when 
backtracking  from  node  3,  we  move  to  the  parent  rode  2,  and  to  the  right 
to  node  4  by  setting  =  0  .  However,  when  backtracking  from  node  6, 

we  move  back  to  node  4,  then  back  to  node  2,  then  back  to  node  1,  and  to 
the  right  to  node  7  by  setting  x^  =  0  . 

On  the  other  hand,  if  the  partial  solution  cannot  be  fathomed, 

we  branch  to  the  left  and  augment  $  by  fixing  a  free  variable  x.,  at 

x,  l  j 

1  (based  on  the  branching  rule),  and  then  we  try  to  fathom  the  resulting 
partial  solution.  In  addition  to  the  one  variable  selected  by  the 
branching  rule,  some  other  free  x. ,  variables  can  also  be  fixed  at 

ij 

0  or  1  according  to  the  application  of  rules  other  than  the  branching 

rule.  Note  that  this  can  also  happen  when  backtracking,  i.e.,  when 

has  been  fathomed  and  we  backtrack  and  move  to  the  right  by  setting 

the  appropriate  x..  variable  to  0  . 
ij 

Let  us  consider  examples  of  both  situations,  i.e.,  when 
has  not  been  fathomed  and  when  has  been  fathomed.  In  Figure  4a 

we  cannot  fathom  (i.e.,  S  at  node  1),  so  we  move  to  node  2  by 


*1 
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augment ing  S,  by  fixing  =  1  based  on  the  branching  rule,  and 

by  fixing  =  0  and  x^  =  0  based  on  the  application  of  the  other 

rules.  Similarly,  we  move  from  node  2  to  node  3  by  augmenting  by 

fixing  x+,  =  l  .  As  an  example  of  backtracking,  wher  we  fathom  , 

we  move  back  lo  the  parent  node  2,  and  to  the  right  to  node  4,  getting 

a  now  partial  solution  S,  by  replacing  x, 0  =  1  with  x, 0  =  0  ,  and 

■a  *4  2.  4  Z 

further  augmenting  it  by  fixing  x(>.  =0  and  x.,  -  1  based  on  the 

^.1  il 

application  of  the  other  rules. 

Computationally,  the  storage  and  update  of  partial  solution 
is  easily  accomplished  by  considering  Figure  4b.  If,  at  a  given  node, 
the  partial  solution  has  not  been  fathomed,  e.g.  ,  at  node  4, 

determine  the  next  branching  variable  by  using  the  branching  rule, 
i.e.,  x.-,^  ,  and  augment  by  adding  204  as  the  last  entry.  Also, 

augment  with  any  other  free  variables,  if  appropriate, 

depending  on  the  application  of  the  other  rules.  Now,  consider  the 
case  where  the  partial  solution  has  been  fathomed,  e.g.,  at  node  6, 

and  we  backtrack;  starting  with  the  last  entry  in  .  we  consider  one 

entry  at  a  Lime,  going  backwards,  until  we  find  a  posi  ive  number  which 
is  not  underlined.  In  our  example,  it  is  103.  In  other  words,  we  must 
branch  to  tb:?  right  by  fixing  X  -  0  ,  i.e.,  we  replace  103  with  -  103 

and  we  are  at  node  7,  Should  we  find  that  we  have  no  positive  number, 
the  procedure  terminates  since  we  are  back  at  the  root  node  and  the 
right,  branch  has  already  been  explored.  This  happens  when  backtracking 
from  node  7. 

In  the  branc h-and-bound  procedure  we  generate  a  sequence  of 
partial  solutions  as  we  move  from  one  node  to  another.  This  sequence  is 
noil-redundant,  in  the  sense  that  no  completion  of  a  partial  solution 
ever  duplicates  a  completion  of  a  previous  partial  solution  that  has 


been  fathomed. 


T 
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Since  one  of  the  x  values,  for  each  j  ,  must  be  1  ,  a  total  of 

ij 

(2m-l)  nodes  are  theoretically  possible  for  complete  enumeration. 

However,  most  of  the  solutions  may  be  infeasible  because  of  the  capacity 
constraints.  The  branch-and-bound  procedure,  thro*.gh  a  judicious  choice 
of  branching  variables,  and  elimination  of  certain  infeasible  and  non- 
optimal  assignments  through  various  rules,  turns  out  to  be  a  practical 
and  computationally  efficent  algorithm.  The  various  components  of  this 
procedure  are  described  next.  Detailed  procedural  steps  and  the  solution 
of  a  test  problem  will  be  covered  in  Chapter  4. 

3 . 1  Bounds 

3.1.1  Lower  Bound 

At  a  given  node  Z  in  the  branch-and-bound  tree,  a  lower  bound 
(LOWS)  is  obtained  by  solving  relaxed  problem  (PR^)  . 

LOWB  -  Z(PR£)  (40) 

Recall  that  problem  (PR£)  is  very  easy  to  solve  by  considering  the 

minimum  over  those  j’s  for  which  x,.  is  not  fixed  at  1  , 

1J  *  1  j 

i.e.,  jeW  p  where  W  is  the  complement  of  W  defined  by  expression  (32). 

Z(PR  )  =  l  c  +  Z_  min  c  +  FC0  ,  (41) 

jeW  jew  i  _  iJx'  ^ 

(i, j)es 


where  c  ^  given  by  expression  (38),  i.e., 


C,-<0  =  a.4  +  “  b,,  O-  -  6,.o)  r<f,,  • 


ij*.  ij 


M'  ijk 


(38) 


and  the  fixed  cost  (FC0)  is  given  by  expression  (39),  i.e.. 


FCt  •  *  V  \  • 

k 


(39) 


where  6^  is  given  by  expression  (34). 


X 


I 


Note  that  it  none  of  the  x.  variables  is  fixed  at  1  ,  as  is 

ij 

generally  the  ease  at  the  root  node,  then  all  6.  ,  =  0  ,  and,  therefore, 

k_L 

P 

FC  -  0  ,  and  c . .  -  a . .  +  L  b  r  .  Z(PR  )  is,  then,  simplv  the 

1  iJ1  iJ  k=1  k  ljk  1  ’  *  ' 

middle  part  ot  expression  (41).  Vve  use  the  term  ’’generally’1  because 
it  is  possible  that  the  capacity  rule  could  force  certain  x 

ij 

variables  to  1  (or  0)  at  the  root  node,  prior  to  solving  the  relaxed 
problem  (PR^)  . 

3.1.2  Upper  Bound 

At  any  given  node  £  ,  let  X  =  {x  }  represent  the  solution  of 
problem  (PR.)  .  Lf  this  solution  is  feasible  for  problem  (P) ,  i.e., 
if  X  satisfies  the  capacity  constraints  (5)  or  (5’) 

•  V  dijk  xij  -  \  yk  Vk  ,  («) 

xijex 

where  y.  =  ]  if  X  d  .  . ,  x  >  0  , 
k  i  j  ljk  lj 

X..CX  (43) 

=  0  otherwise  , 


then  the  value  of  problem  (P)  corresponding  to  this  solution 
gives  an  upper  bound  (UPB) : 


XL  a  ,  .  x .  .  +  L  b  v  , 
.  .  lj  ij  .  k  k 

i  J  J  J  k 

x .  .ex 

-u 


where  y  is  defined  by  (43)  . 


3.1.3  Test  Upper  Bourd 


A  current  Lowest  upper  bound  is  retained  as  the  lest  upper  bound 
(BUB),  the  correspond  trig  solution  X  representing  the  incumbent  solution. 
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The  branch-and-bound  procedure  is  initiated  by  assuming  a  very  large 
value  as  the  best  upper  bound,  and  is  replaced  by  better  (lower)  values 
as  the  procedure  continues. 

A  positive  fractional  value  c  can  be  specified  if  a  sub-optimal 
solution  is  acceptable.  For  example,  for  £  =  0.001  ,  the  resulting 
solution  value  is  guaranteed  to  be  within  0.1  percent  of  the  optimal 
solution  value.  When  £  is  non-zero,  the  adjusted  best  upper  bound  (BUBS) 
is  defined  as: 

BUBS  =  BUB/ (1  +  £)  .  (45) 

Obviously  when  0=0,  BUBS  =  BUB  . 


3. 2  Facility  Usage  Rule 

This  rule  is  used  to  identify  facilities  forced  into  usage  at  a 
given  node  £  and  hence  fix  corresponding  free  variables  y^_  at  1  . 


For  a  partial  solution  S.  ,  define 


V  ■  dijk  lf  j£"  ' 


-  min  d . . .  if  j  eW  . 
±  Uk 

(i, j)e§ 


(46) 


The  facility  usage  rule  states  that  for  any  facility  k  ,  where 
y^  is  not  already  fixed  at  1  ,  if  E  d,^  >  0  ,  tnen  facility  k  is 

j  J 

forced  into  usage  and,  therefore,  should  be  fixed  at  1  . 

This  rule  is  applied  at  every  node  prior  to  applying  the  capacity 
rule.  In  other  words,  this  rule  is  applicable  to  capacitated  as  well  as 
uncapacitated  problems. 


3 . 3  Capacity  Rule 

This  rule  is  designed  to  "exclude”  infeasible  assignments  prior 
to  solving  the  relaxed  problem  (PR^)  •  This  is  done  by  exploiting  the 
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relationship  between  the  capacities  requited  (d  .  .  and  the  capacities 

1  j  k 


available  (su)  lor  a  given  partial  solution  of  problem  (P), 


The  capacity  rule  states  that  for  a  facility  k  and  an  activity 


j  , "exclude" a  free  x..  variable  (i.e.,  fix  it  at  0)  for  which 

ij 


(d  .  . .  -  d.L  )  >  (s.  -  E  d  0 ) 
ijk  jk  i  k  .  jkx 


(  i  *  j  )  ‘  S 


(47) 


where  d . ,  ,  is  defined  by  expression  (46).  The  right-hand  side  of  this 
j  kx 


inequality  (47),  when  positive,  represents  the  available  capacity  at 
facility  k  .  The  left-hand  side  shows,  for  a  given  j  ,  the  difference 


between  a  d...  cor respond ing  to  a  free  x..  variable  and  d..«  . 

ijk  r  a  ij  jkx 


If,  for  a  specific  d  . ,  ,  this  difference  is  more  than  the  available 

ijk 


capacity,  the  corresponding  free  x..  variable,  if  fixed  at  1  ,  would 


result  in  an  infeasible  solution.  Thus,  by  looking  ahead,  we  can  exclude 


such  a  free  variable  by  assigning  it  a  value  of  0  . 


Note  that  if  the  right-hand  side  of  expression  (47)  is  negative, 
then  any  completion  of  such  a  partial  solution  will  be  infeasible  and 
we  backtrack  in  our  branch-and-bound  procedure. 


The  capacity  rule  is  applied  to  all  the  facilities  by  considering 
one  facility  at  a  time.  The  cycle  of  examining  all  the  facilities 
continues  until  no  more  assignments  can  be  excluded.  During  the  course 


of  application  of  this  rule,  if  all  but  one  of  the  free  x^  variables 


have  been  excluded  (fixed  at  0)  for  a  given  j  ,  then  that  particular 


x..  variable  is  fixed  at  1  because  of  constraints  (2),  i.e.,  each 
ij 


activity  j  must  be  assigned  to  one  and  only  one  design  i  .  The 


partial  solution  is  updated  accordingly  to  reflect  the  x  variables 


fixed  at  0  or  1  due  to  the  application  of  the  capacity  rule. 


The  capacity  constraints  for  an  uncapacitated  problem  are  not 

Hence,  the  capacity  rule  is  useful  only  for  capacitated  problems 


ac t ive. 
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3 . 4  Branch ing  Rule 

This  rule  provides  the  choice  of  the  x, .  variables  on  which  to 

ij 

branch.  If  the  partial  solution  at  a  given  node  £  is  not  fathomed,  we 
branch  further  by  fixing  a  free  variable  at  1  and  moving  to  the 

left  branch  node. 


According  to  the  branching  rule  the  choice  of  the  branching 

variable  depends  on  the  c  .  values  and  is  such  that  the  corresponding 

ij  * 

x_  ,  if  perturbed,  has  the  maximum  impact  on  the  optimal  value  of 
problem  (PR^)  • 

For  a  given  j  ,  define  c  Q  ,  the  minimum  permissible  c  i 

ljx, 


and  c,  , . 

,  the  second  smallest  permissible 

Cii£  ’  i,e” 

ci1j^  =  mjn  cij£ 

for 

jeW  and 

(i,j)es 

(48) 

and 

c±2j£  =  min  Cij£ 

for 

j  eW  and 

(i, j)e§ 

(49) 

For 

each  jeW  ,  define 

V 

=  Ci2j&  " 

Ci13^  ‘ 

(50) 

Our  branching  rule  states  that  a  free  x^  variable  corresponding  to 

c.  . Q  such  that  D,0  is  maximized  over  all  j  ,  is  selected  as  the 
i  ^  3  ^  J  ~ 

next  branching  variable  and  assigned  a  value  of  1. 


3 . 5  Bounding  Rule 

This  rule  is  designed  to  "exclude’*  certain  non-optimal  assignments . 
These  assignments  cannot  lead  to  an  optimal  solution  as  we  branch  from 
one  node  to  the  next  left  branch  node. 


The  bounding  rule  states  that  a  free 


variable  should  be 


excluded  (by  assigning  it  the  value  0)  for  which 


and 


(c..,  -  c.  .  ,>)  >  (BUBS  -  LOWB)  for  jeW 


alj  ' 


where  c  ,  ,  BUBS,  and  LOWB  are  given  by 
(40),  respectively. 


and  (i,j)cS 
expressions  (48) , 


(51) 

(45), 


Thus,  by  looking  ahead,  we  exclude  those  assignments  which  will 
provide  lower  bounds  higher  than  BUBS. 

The  bounding  rule  is  applied  to  each  jeW  just  prior  to  selecting 

the  x  variable  for  branching  to  the  left, 
ij 

As  in  the  case  of  the  capacity  rule,  if  the  bounding  rule  results 
in  excluding  (fixing  at  0)  all  but  one  of  the  free  x  variables  for 

a  given  jCW  ,  then  that  particular  x..  variable  is  fixed  at  1  .  Also 

the  partial  solution  is  updated  accordingly  to  reflect  the  x  variables 

fixed  at  0  or  1  due  to  the  application  of  the  bounding  rule. 


3 . 6  Bac kt racking  Rules 


If  a  partial  solution  at  a  given  node  has  beer,  fathomed,  we  back¬ 
track.  The  backtracking  rules  are  typical  of  a  branch-nnd-bound  procedure. 
In  addition,  the  application  of  the  capacity  rule  and  the  bounding  rule 
can  lead  to  backtracking.  The  criteria  for  backtracking  include  the 
fu 1  low  ing. 


(a)  When  applying  the 
capacity  given  by 
(47)  is  negative, 

then  backtrack. 


capacity  rule, 
the  right-hand 


l.  e. 


-  I 
j 


if  the  available 
side  of  inequality 


(b)  If  LOWB  "■  BUBS,  then  backtrack.  Otherwise  compute 
uru  if  the  solution  is  feasible  in  problem  (P)  . 

Then  update  BUB  and  BUBS  if  UPB  <  BUB;  and  back¬ 
track  if  LOWB  =  BUBS. 

(c)  If  further  branching  is  not  possible,  then  backtrack. 

This  can  happen  due  to  the  capacity  rule,  the  bounding 
rule,  or  the  branching  rule  if  the  updated  partial 
solution  is  such  that  no  further  branching  is  possible, 
i.e.,  x,.  variables  are  fixed  at  J  for  all  j  ,  or 

equ ivu lent ly ,  W  *  $ 


IT 
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When  any  of  the  backtracking  criteria  apply,  we  backtrack  to 
the  parent  node  and  move  to  the  right  branch  node  (if  the  right  branch 
has  not  already  been  explored)  by  fixing  the  appropriate  x^  variable 

at  0  .  If  the  right  branch  has  already  been  explored,  we  continue  back¬ 
tracking  to  a  parent  node  where  we  can  move  to  a  right  branch  node. 

The  branch-and-bound  procedure  terminates  when  we  backtrack  to  the 
root  node  and  find  that  the  right  branch  node  has  already  been  explored. 


4.  COMPUTATIONAL  STEPS  AND  THE  COMPUTER  PROGRAM 


A  computer  program  called  ZIPCAF  (an  acronym  for  Zero-one  Integer 
Program  for  mul t iac liv  ity  multifacility  Capacity-.  onstrained  Assignment 
Problems)  implementing  the  branch-and-bound  methouology  has  been  developed. 

Detailed  procedural  steps  and  guidelines  to  use  the  computer 
program  are  described  in  a  separate  document  [Chhabra  and  Soland  (1980)] 
titled  "Program  Description  and  User's  Guide  for  ZIPCAP — a  Zero-one 
Integer  Program  to  solve  multiactivity  multifacility  Capac itv-constrained 
Assignment  Problems."  Specifically,  the  document  includes: 

.  Problem  formulation  (P)  and  potential  areas  of 
appl ica t ion 

.  Overall  flow  diagram  and  detailed  procedural  steps 
for  the  computer  program 

.  Program  listing  and  dictionary  of  the  symbolic  names. 

The  listing  includes  extensive  use  of  comment  cards 
to  explain  various  computational  steps. 

.  User  information  including 

schematic  diagram  of  the  deck  structure, 

-  detailed  instructions  for  the  job  control  (JCL) 
cards,  program  parameter  card,  program  options 
card,  and  the  various  other  input  data  cards. 

.  Three  test  problems  to  demonstrate  the  use  of  the 
program.  The  dispJay  Includes  coded  input  and 
annotated  outputs  reflecting  the  use  of  selected 
program  options. 

As  mentioned  ear  ier,  ZIPCAP  is  primarily  designed  for  capacitated 
problems.  However,  uncapacitated  problems  can  be  solved  as  a  special 
case,  and  this  is  demonstrated  by  including  an  uncapacitated  test 
problem. 
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Because  of  the  extensive  coverage  of  the  program  description  and 
user  guidelines  in  the  above  document,  this  chapter  provides  only  an  over¬ 
view  of  the  computer  program,  including  an  overall  flow  diagram,  and  a 
summary  of  the  program  options,  in  order  to  provide  continuity  in  this 
document.  In  addition,  a  step-by-step  description  Ov  a  test  problem  is 
presented  to  demonstrate  the  use  of  the  various  components  of  the  branch- 
and-bound  methodology.  The  computer  printout  showing  step-by-step  details 
is  obtained  by  use  of  one  of  the  program  options.  The  use  of  this  option 
to  display  detailed  steps  in  this  document,  in  fact,  complements  the  use 
of  the  various  options  demonstrated  in  the  other  document. 

4 . 1  The  Program 

Figure  5  presents  a  simplified  flow  diagram  of  the  branch-and- 
bound  procedure.  The  major  computational  steps  for  the  computer  program 
are  numbered  in  circles.  These  steps  are  essentially  based  on  the 
methodology  components  described  in  the  previous  chapter.  A  step-by-step 
description  has  been  included  in  the  other  document  [Chhabra  and  Soland 
(1980) ] . 

The  computer  program  ZIPCAP  is  written  in  FORTRAN  IV,  and  has 
been  developed  and  tested  on  an  IBM  3031  at  the  George  Washington 
University.  The  program,  comprising  about  480  lines  is  currently 
dimensioned  for  a  maximum  problem  size  of  35  designs  (m^ ,  35  activities 
(n)  and  30  facilities  (p) .  The  program  size  to  execute  a  problem  has 
two  components:  one,  due  to  the  program  itself,  comprising  173  K  bytes, 
and  the  other  dependent  on  the  dimensions  of  the  arrays  given  by  the 
following  functional  relationship. 

f(m,n,p)  =  4[(p+4)mn  +  (m+5)p+9n]  bytes 

The  computer  program  listed  in  the  other  document  has  since  been 
further  improved.  The  basic  improvement  has  been  the  addition  of  the 
facility  usage  rule.  This  rule,  as  described  in  Chapter  3,  is  applied 
both  to  capacitated  and  uncapacitated  problems  just  before  the  applica¬ 
tion  of  the  capacity  rule.  For  completeness  of  this  document,  a  revised 
program  listing  is  included  in  Appendix  A.  It  may  te  mentioned  that 
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Figure  5. 

Simplified  flow  diagram  for  the  branch-and -bound  procedure 
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the  revised  program  solves  the  test  problems  included  in  the  other  document 
more  efficiently  —  in  less  time  and  in  fewer  noc  es  (with  an  average 
reduction  in  nodes  of  31  percent) .  The  improvement  in  efficiency  seems 
to  result  from  the  "multiplicative”  effect  of  the  various  rules.  Another 
improvement  made  is  that  the  computer  printout  always  displays  the  node 
number  (IBNOD)  at  which  the  best  upper  bound  changes  (improves)  and  the 
corresponding  values  of  the  best  upper  bound  (BUB)  and  the  adjusted  best 
upper  bound  (BIBS) . 

ZIPCAP  provides  numerous  options  to  the  program  user.  These  options, 
described  in  the  other  document,  are  summarized  in  Table  3. 

Option  ICAPR,  the  capacity  rule,  is  automatically  skipped  by  the 
program  when  solving  an  uncapacitated  problem.  Option  ISTEP,  the 
intermediate  steps'  listing,  even  when  skipped,  provides  information  on 
the  total  number  of  nodes  explored.  A  summary  listing  provides  necessary 
information  tc  construct  the  branch-and-bound  tree,  whereas  a  detailed 
listing  of  the  intermediate  steps  is  useful  when  changing  or  debugging 
the  program. 

Option  EPS,  the  optimal/suboptimai  solution,  provides  the 
flexibility  of  obtaining  a  suboptimal  value  guaranteed  to  be  within 
a  specified  fraction  of  the  optimal  value.  The  resulting  solution  may 
he  suboptimal  but  could  provide  a  considerable  saving  in  terms  of 
exploring  fewer  nodes  in  comparison  to  those  necessary  for  obtaining 
an  optimal  solution. 

Option  ET,  by  providing  important  information  at  a  specified 
elapsed  time,  is  useful  in  a  situation  where  the  total  time  allocated 
to  solve  a  problem  may  not  be  sufficient  and  the  program  terminates 
before  verifying  an  optimal  solution.  The  information  provided  by 
this  option  includes  an  updated  partial  solution  showing  the 

variables  fixed  at  0  or  1,  at  the  current  node  being  explored  at 
the  specified  tire  ET.  By  looking  at  the  first  few  variables  displayed 
in  the  partial  solution  of  the  current  node,  it  is  possible  to  assess 
the  extent  ot  the  branch-and-bound  tree  explored  until  time  ET.  For 


SUMMARY  OF  ZIPCAP  OPTIONS 


Skip  this  option 


example,  in  view  of  the  terminology  in  Figure  4b  (Chapter  3),  if,  at  an 

arbitrary  node,  the  first  term  of  the  partial  solution  is  positive, 

I.e.,  the  x  .  variable  has  value  1,  then  we  are  still  in  the  left 
iJ 

half  of  the  total  branch-and-bound  tree.  If  the  first  term  is  negative, 
i.e.,  the  x^.  variable  has  value  0,  then  we  ar^.  in  the  right  half  of 

the  total  branch-and-bound  tree  and  have  explored  half  of  the  total 
(theoretical)  solutions  corresponding  to  the  left  half  of  the  tree. 

If  the  first  two  terns  are  negative,  i.e.,  the  fi-st  two 

variables  have  value  0,  then  one  quarter  of  the  total  (theoretical) 
solutions  remain  to  be  explored,  since  we  are  in  the  next  right  half  of 
the  right  half  of  the  total  branch-and-bound  tree,  -.s  illustrated  in 
Figure  6. 


Left  half  of 
the  "total" 
branch-and- 
bound  tree 


|  Each  t egment  represents  1/4 

l  of  the  total  (theoretical) 

solutions 

Each  segment  represents  1/2 
Right  half  of  of  total  (theoretical) 

the  "total'*  solutions 

branch-and- 
bound  tree 


Figure  6 

Illustration  for  estimating  the  extent  of  the  branch 
and-bound  tree  explored 


Recall  f rou’  Chapter  3,  that  a  total  of  (2m-l)n  nodet  are  theoretically 
possible.  Thus,  it  the  first  g  [g  <  (m-l)n]  terms  at  an  arbitrary 

node  are  negative,  then  theoretically  about  [  (2n— l)n/^g]  nodes  remain 
to  be  explored. 


4 . 2  An  Illustrative  Example 

We  consider  a  capacitated  test  problem  with  five  designs  (m)  , 
four  activities  (n) ,  and  eight  facilities  (p)  to  demonstrate  the  use 
of  the  branch-and-bound  procedure  and  the  computer  program. 

The  computer  printout  for  this  problem  shoving  step-by-step 
details  for  a  couple  of  nodes  is  presented  in  Appendix  B. 

As  shown  in  the  beginning  of  the  printout,  the  options  selected 

are: 


.  IINPT  =  1,  i.e.,  list  the  input  data 

.  ICAPR  =  1,  i.e.,  use  the  capacity  rule 

.  ISTEP  =  2,  i.e.,  list  detailed  intermediate  steps 

.  IUNCAP  =  1,  since  this  is  a  capacitated  problem 

EPS  =0.0  implying  that  an  optimal  solution  is  desired 

.  ET  =  0.0  since  a  detailed  listing  of  intermediate  steps 

will  be  available. 


Following  the  listing  of  the  options,  input  data  listed  for  the  problem 
include  variable  costs  a..  ,  fixed  costs  b,  ,  available  capacities 

,  and  capacities  required  d  _ ^  .  The  values  are  generated 

by  the  computer  program. 

The  computer  program  follows  the  procedural  steps  marked  in  the 
flow  diagram  presented  in  Figure  5.  These  steps,  llong  with  the 
relevant  terminology  used  in  the  computer  printout,  are  described  below 
for  a  couple  of  nodes,  followed  by  a  complete  branch-ana -bound  tree  for 
this  problem.  As  mentioned  earlier,  a  dictionary  of  the  symbolic  names 
used  in  the  computer  program  is  included  in  the  other  document. 
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Node  1 

Step  J:  Initialize. 

Initialize  BUB  =  9999999.0,  and  since  EPS  =  0.0,  BUBS  =  BUB. 

Also  S  =  0  and  W  =  <f>  .  In  the  computer  printout,  vector 

FIX(J)  represents  the  s£l  W,  and  matrix  CX(I,J)  ;epresents  both,  fixed 

and  free  x  .  variables.  In  the  CX(I,J)  matrix,  an  x.,  variable  fixed 
J-J  iJ 

at  1  or  0  is  represented  as  1  or  2,  respectively,  and  a  free  x  variable 

is  represented  by  the  value  0.  Initially,  all  the  x_  variables  are  free 

as  shown  by  matrix  CX(I,J)  in  the  printout. 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l, 2 , . * . , 8  . 

In  the  printout,  MIND(J)  represents  d  ^  defined  by  expression  (46), 
and  MINSD  represents  I  d.^  .  As  shown  in  the  printout,  MINSD  is  0  for 

.i  3 

k=l,2,...,8  ,  and  so  the  facility  usage  rule  does  not  orce  any  facilities 

into  usage;  and  as  shown  by  matrix  CX(I,J)  for  k=l,2,...8  ,  the  capacity 

rule  does  not  fj.x  any  x,.  variables. 

J-J 

Step  3:  Solve  the  relaxed  problem  (PR^)  * 

In  the  printout  FLB(K)  represents  6^  *  given  by  expression  (34), 

lor  computing  FCh  ,  and  C(I,J)  represents  c . . Q  defined  by  expression 

ij  ~ 

(3B).  Being  at  the  root  node,  %  =  1.  Further  the  solution  of  problem 
(PR.)  ,  i.e.,  X  =  {x..}  is  shown  in  the  printout  by  SOLX(J)  which 

for  (PR.)  is  X  =  { x ,  .  »  x/0  =  x0.5  =  x, ,  =1}  . 

1  41  42  23  4*4 

Step  4:  Compute  the  lower  bound. 

The  expressions  (40)  and  (41),  i.e., 

LOWB  *  Z(PR^)  (A0) 

-  X  c..n  +  X  min  c  +  FC  . 

.  fT  xj £  .  .  ij£  *  (A  1 ) 

jcW  jr.W  i 

(i, j)eS 


55 


are  represented  in  the  printout  as 

LOWB  -  M1NSC  +  FC 

=  729839.3125  +  0  =  729839.3125 


Step  5:  Compare  LOWB  with  BUBS. 
Since  LOWB  <  BUBS,  go  to  Step  6 


Step  6:  Check  if 
expression  (42)  is 

Z  Z 

1  j 


solution  X  is  feasible 
satisfied . 


dijk  xij 


-  Vk 


X.  .EX 
1J 


in  problem 

Vk 


(P), 


i.e.  , 


(42) 


In  the  printout,  NSUMD  represents  the  left-hand  side  of  this  inequality, 
and  for  each  k  ,  the  capacity  constraints  are  satisfied. 

Step  7:  Compute  the  upper  bound. 

UPB  is  given  by  expression  (44),  i.e., 


UPB  =  Z  Z  a.  .  x.  .  +  Z  b.  y.  (44) 

ij  1J  1J  k  kk 

X..CX 

U 


In  the  printout,  the  corresponding  expression  is  represented  as 

UPB  =  NS UMA  +  FCUB 

=  678,502  +  101,000  =  779502.0  . 

Step  8:  Compare  UPB  with  BUB. 

Since  UPB  <  BUB,  go  to  Step  9. 

Step  9:  Set  BUB  =  779502.0  .  Since  EPS  -  0.0,  BUBS  -  BUB. 

Since  LOWB  <  BUBS,  go  to  Step  10. 

Step  10:  Left  branching  is  possible  since  W  -  <P  as  shoxm  by 
vector  FIX(J);  go  to  Step  11. 

Step  11:  Apply  the  bounding  rule  and  the  branching  rule. 


ij 


According  to  our  bounding  rule,  a  free  x 


variable  is  excluded 


t 
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(fixed  at  0)  for  which 

(c. -  c.  .f)  >  (BUBS  -  LOWB)  for  jeW  and  (i,j)eS  (51) 

ijX/  l^x 

For  x  ,  (210,381.4375-145,201.5)  >  (779 , 502 . 0-7 29 , 839 . 3125) . 

This  also  holds  for  x„^  and  x,  ,  ,  i.e.,  the  bounding  rule  results 
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in  fixing  x^  ,  x^  ,  and  x^  at  0.  This  is  shown  in  the  printout 

by  matrix  CX(I,J)  where  the  corresponding  variables  have  been  assigned 
the  value  2  because  of  the  bounding  rule. 

The  branching  rule  directs  us  to  select  a  free  x..  variable 
°  1J 

corresponding  to  c.  .  0  for  which  D.c  =  c.  . «  -  c.  is  maximized 

J  **■  1^3^ 

over  all  j  .  In  the  printout,  c.  . Q  ,  c .  .Q9  and  D . Q  are  represented 

by  NMINC(J),  MINC(J)  and  DIFBR(J),  respectively.  Since  ^  is  the  maxi¬ 
mum,  x.,,  is  selected  as  the  next  left  branching  variable.  This  is  shown 
42 

in  t fie  printout  by  BR1  and  is  represented  as  (100  i  i-  .; )  e.g.,  402  . 


Using  the  terminology  employed  in  Figures  4a  and  4b,  the 

variables  fixed  at  0  or  1  in  the  partial  solution  will  be  shown 

as  =  {-  103,  -  303,  -  104,  402}  .  In  the  computer  printout,  vector 

STX  displays  the  x_  variables  fixed  at  0  or  1 .  The  representation 

of  the  variables  is,  however,  somewhat  different.  An  x..  variable 

fixed  at  0,  due  to  any  rule,  is  shown  as  -  (100  i  +  j)  -  1,000,000, 
e.g.  ,  x,~  is  shown  as  -  1,000,103;  an  x..  variable  fixed  at  1  due 

to  the  branching  rule  is  represented  as  (100  i  +  j),  e.g.,  x^  as  402; 

and  an  x..  variable  fixed  at  1  due  to  a  rule  other  than  the  branching 

rule  is  shown  as  (100  i  +  j)  +  1,000,000,  e.g.,  x  is  represented  as 


1,000,203. 

In  the  printout,  vector  STX  represents  updated  partial  solution 


We  now  move  Co  Node  2. 


Node  2 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  are  displayed  in  the 
pr intou  t . 


Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l, 2, ... ,8  . 


As  shown  in  the  printout,  MINSD  (representing  £  d  )  ,  being 

j  3 

positive  for  1=1, 2, 3, 4,  and  5  ,  these  facilities  are  forced  into 
usage.  Further,  for  k=4  ,  expression  (47)  holds  for  and  , 


i.e., 


(180-0)  >  (200-30) ,  and 

(180-0)  >  (200-30),  respectively. 

As  shown  by  matrix  CX(I,J)  in  the  printout,  these  two  variables  are 
excluded  (fixed  at  0)  by  the  capacity  rule.  Since  the  capacity  rule 
results  in  fixing  at  least  one  variable  in  the  first  cycle,  another 
cycle  is  repeated  as  displayed  in  the  printout.  The  second  cycle 
does  not  fix  any  more  variables.  Vector  STX  is  updated  accordingly. 

Step  3:  Solve  the  relaxed  problem  (PR^)  * 

represented  by  FLB(K)  ,  c^j2  represented  by  matrix  C(I,J), 
and  solution  y  represented  by  SOLX(J)  are  displayed  in  the  printout. 
Step  4:  Compute  LOWB. 

LOWB,  from  the  printout,  is  equal  to  749011  4375. 


Step  5:  Compare  LOWB  with  BUBS. 

Since  LOWB  <  BUBS,  go  to  Step  6. 

Step  6:  Check  if  solution  X  is  feasible  in  (P) . 

In  the  printout,  for  k=4  ,  NSUMD  =  290  >  200  ,  i.e.,  expression 
(42)  is  not  satisfied,  and  we  go  to  Step  10. 

Step  10:  As  shown  by  vector  FIX(J) ,  left  branching  is  possible 
and  we  go  to  Step  11. 


58 


Step  11:  Apply  the  bounding  rule  and  the  branching  rule. 

As  displayed  by  matrix  CX(I,J)  in  the  printout,  the  bounding 
rule  results  In  fixing  and  at  0.  Now,  for  j=4  ,  except 

for  ,  all  the  x„  variables  are  fixed  at  0;  therefore  , 

is  fixed  at  1.  This  is  reflected  by  matrix  CX(I.J),  and  vector  FIX(J). 
Vector  STX  is  updated  accordingly. 

The  branching  rule  selects  x^  as  the  next  branching  variable. 

This  is  shown  in  the  printout  by  BR1,  and  vector  STX  is  updated 
accordingly. 

We  now  iove  to  Node  3. 


Node  3: 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  are  displayed  in  the 
printout . 

Step  2;  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k=l, 2, ... ,8  . 


The  facility  usage  rule  forces  facilities  1  to  5,  and  8  into 
usage.  For  k=4  ,  the  capacity  rule  excludes  and  x«^  ,  i.e., 

fixes  them  at  0;  and  for  j=3  ,  all  but  being  fixed  at  0, 

is  fixed  at  1.  This  is  displayed  in  the  printout  by  matrix  CX(I,J) 
and  vector  FIX(J).  Vector  STX  is  updated  accordingly. 


Although  the  capacity  rule  has  fixed  at  least  ona  x^  variable 

during  the  initial  cycle,  another  cycle  is  not  necessary,  as  displayed 
by  vector  FIX(J)  which  represents  set  W,  since  we  have  an  x^  variable 

fixed  at  1  for  each  of  the  n  columns  (activities) . 


Step  3:  Solve  the  relaxed  problem  (PR^)  • 

SOLX(J)  displays  the  solution  for  the  relaxed  problem. 


Step  4:  Compute  LOWB. 


LOWB,  shown  in  the  printout,  is  equal  to  7 7 9502. 0  . 
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Step  5:  Compare  LOWB  with  BUBS. 

Since  LOWB  =  BUBS,  go  to  Step  12. 

Step  12:  Backtrack. 

We  backtrack  by  moving  to  the  parent  Node  2,  ard  branching  to  the 
right  by  setting  =  0  (since  the  right  branch  has  not  yet  been  explored). 

In  the  printout,  this  is  accomplished  by  observing  the  last  entry  in  vector 
SIX,  and  moving  backwards,  one  entry  at  a  time,  until  we  find  a  positive 
entry  without  1,000,000  added  to  it.  The  corresponds ng  variable 

is  fixed  at  0,  and  we  move  to  the  right  branch  node.  Matrix  CX(I,J),  vector 
FIX(J)  and  vector  STX  are  updated  accordingly.  As  displayed  in  the  print¬ 
out,  entry  401  in  vector  STX  is  such  an  entry,  and  variable  is  fixed 

at  0  for  branching  to  the  right.  This  is  shown  in  the  printout  by  BR0  as  401. 
The  updated  vector  STX  is  also  displayed. 

We  now  move  to  Node  4. 


Node  4 

The  updated  matrix  CX(I,J)  and  vector  FIX(J)  are  displayed  in  the 
pr intout . 

Step  2:  Apply  the  facility  usage  rule  and  the  capacity  rule  for 
k*l , 2 , . . . , b  . 


As  displayed  in  the  pr intou t ,  f or  k=4,  MINSD=230  >  200  ,  i.e., 
the  right-hand  side  of  inequality  (47),  (s^  -  E  cL^)  <  ^  ’  anc* 

according  to  our  backtracking  rules,  we  backtrack,  i.e.,  go  to  Step  12, 


Step  12:  Backtrack. 


We  backtrack  to  the  parent  Node  2,  and  since  the  right-hand  branch 
has  already  been  explored,  backtrack  to  Node  1  and  to  the  right-hand 
branch  by  fixing  x ^  to  0  •  This  is  shown  in  the  printout  by  BR0  as  402, 

and  vector  STX  is  updated  accordingly. 

We  now  move  to  the  next  node,  i.e..  Node  5. 
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Branch-and-Bound  Tree 

We  continue  the  branch-and -bound  procedure  from  one  node  to 
another  until  we  backtrack  to  the  root  node  and  find  that  the  right 
branch  has  already  been  explored.  The  procedure,  then,  terminates  and 
the  solution  corresponding  to  the  best  upper  bound  is  the  optimal  solu¬ 
tion. 

For  this  problem,  a  total  of  nine  nodes  are  explored  and  the 
optimal  value  equals  779502.0.  The  optimal  solution  is  x^=  x^2~ 

x23  =  X44  =  1  and  yl  =  y2  =  y3  =  y4  =  y5  =  y8  "  1  *  ThiS  iS 
displayed  in  the  computer  printout  on  the  last  page  of  Appendix  B. 

Figure  7a  presents  the  branch-and-bound  tree  fc  r  this  problem, 
and  shows  the  node  numbers,  the  bounds,  and  the  branching  variables. 

In  ordei  to  demonstrate  the  role  of  the  capacity  rule  and  the 
bounding  rule,  Figure  7b  displays  the  variables  fixed  as  0  or  1 

by  these  rules  for  this  test  problem. 

The  cumulative  effect  of  the  various  rules,  including  the  facility 
usage  rule,  the  capacity  rule,  and  the  bounding  rule,  , lakes  the  branch- 
and-bound  procedure  quite  efficient.  Further,  the  storage  and  updating 
of  the  x  variables  fixed  at  0  or  1  is  done  in  a  manner  that  makes 

ij 

utmost  use  of  tne  relevant  information  at  the  preceding  node. 


■aatfMMiiliflii 


LOWB  =  729839.3 
JJPB  =  779502.0  =  BUB 


x ,  ,-,-1 

%/ 

LOWB  =  749013.412) 


X41=V 


779502.0  (3. 


^42=0 


X23  li 


5  )  LOWB  =  747995.3 


6)  LOWB  =  771864. 9 
h-\UPB  =  857825.0 


9)  LOWB  =  793044.0 


LOWB  =  807947.0  vy  cy 

Figure  7a 

Branch-and-bound  tree  for  a  test  problem 


(Test  Problem  with  m=5 ,  n=4 ,  and  p=8) 

Node 

Capacity  Rule 

Bounding  Rule 

1 

■  ■ 

X13=°  *  X33=°  ’  X14=° 

2 

X34=°  *  X54=0 

x2r°  ’  x24=0  ’  x44=1 

3 

x^3~“0  ,  x^-0  ,  X2-j— 1 

5 

X43=°  ’  X34=° 

6 

XH=0  ,  x7^=0  .  x  =0  »  x22=°  ’  X54=t> 

8 

rH 

■1 

*  1 

1 

9 

x53  1 >  x44“0,  x54  x24~* 

Figure  7b 

Variables  fixed  by  the  capacity  rule  and  the  bounding  rule 


5.  COMPUTATIONAL  TEST  RESULTS 


The  computer  program  ZIPCAP  has  been  tested  on  several  problems. 
Although  primarily  designed  for  capacitated  problemr.  (i.e.,  where  the 
capacity  constraints  are  active),  the  program  can  also  be  used  for  solving 
uncapacitated  problems  as  a  special  case.  Since  the  dita  available  for 
capacitated  problems  were  limited,  some  uncapacitated  problems  were  also 
considered  foj  testing  the  program.  (Most  of  the  data  were  furnished  by 
Professor  Pinkur  and  are  related  to  his  work  on  multi-echelon  inventory 
systems. ) 

Table  4  presents  the  test  results  of  ZIPCAP.  In  order  to  verify 
the  optimal  solutions,  the  test  problems  were  also  solved  by  using  the 
0-1  integer  programming  code  RIP30C  [Geoff rion  and  Nelson  (1968)]. 

In  the  table,  the  problem  size  shows  the  number  of  designs  (ra)  , 
activities  (n) ,  and  facilities  (p) .  This  is  equivalent  to  solving  a 
problem  having  ran+p  variables  and  n+p  constraints.  The  elapsed  time 
represents  the  time  in  seconds  to  solve  the  problem,  excluding  the  time 
to  read  and  write  the  input  and  to  write  the  output.  The  total  number 
of  nodes  explored  by  ZIPCAP  for  a  specified  set  of  options  is  also  shown. 

Both  RIP30C  and  ZIPCAP  were  run  on  an  IBM  3031  al  The  George 
Washington  University.  The  last  problem  in  the  table  was  not  run  using 
RIP30C  because  of  the  code’s  capacity  limitation  to  $0  variables  and 
50  constraints. 

The  test  problem  with  m^  ,  n=4  ,  and  p=5  has  three  variations, 
using  different  values  for  the  facility  capacities.  The  data  for  the 
variable  costs  ,  fixed  cost  ,  and  the  capacity  requirements 

d...  are  given  in  the  other  document,  i.e.,  Chhabra  and  Soland  (1980). 
lj  k 

For  the  test  problem  with  m*5  ,  n=4  ,  and  p«=8  .  runs  4a,  4b, 
and  4c  are  the  same  except  for  nc  different  intermediate  steps’  option 
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RIP30C  ZIPCAP 
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(ISTEP)  and  this  results  in  slight  differences  in  the  time  taken  to  solve 
the  problem.  Runs  4a  and  4d  differ  in  that  4d  does  n:>t  use  the  capacity 
rule;  the  resulting  difference  in  the  total  number  of  nodes  explored 
to  reach  the  optimal  value  points  to  the  effectiveness  of  the  capacity 
rule  in  conjunction  with  the  bounding  rule. 

Run  5  shows  the  results  for  an  uncapacitated  problem  with 
m=10  ,  n=8  ,  and  p=8  .  Option  ICAPR  is  not  used  since  the  capacity 
rule  is  not  useful  for  an  uncapacitated  problem. 

Another  uncapacitated  problem  with  m=10  ,  n-30  ,  and  p-8  is 
solved  in  runs  6a  and  6b.  In  run  6a,  the  epsilon  value  (EPS)  is 
specified  as  0.002.  The  solution  value  found  by  exploring  125  nodes 
may  be  subopt Lmal  but  is  guaranteed  to  be  within  +0.2  percent  of  the 
optimal  solution  value.  Run  6b  is  made  with  an  epsilon  value  (EPS) 
of  0.0,  and  the  optimal  solution  value  is  found  in  277  nodes.  A  comparison 
of  runs  6a  and  6b  shows  that  the  number  of  nodes  :.s  less  than  half  for  a 
solution  value  that  may  be  suboptimal  but  very  close  to  the  optimal  solu¬ 
tion  value  as  compared  to  the  number  of  nodes  for  an  optimal  solution 
value. 

In  general,  a  small  difference  between  a  solution  value  that  may 
be  suboptimal  and  the  optimal  solution  value,  translates  into  a 
significant  difference  in  the  corresponding  number  of  nodes  and  the 
solution  time  required. 


6.  FURTHER  CONSIDERATIONS 


It  was  mentioned  in  Chapter  2  that  it  is  possible  to  consider 
alternative  formulations  of  problem  (P) ,  and  also  to  c  msider  choices  of 
Lagrange  multipliers  other  than  v  =  b  with  the  purpose  of  obtaining 

K  K 

“tighter"  bounds  which,  in  turn,  would  further  improve  the  efficiency 
of  the  branch~and-bound  procedure.  These  aspects  will  be  discussed  in 
this  Chapter. 


6.1  Alternative  Formulations 

Problem  (P)  can  be  reformulated  by  adding  additional  constraints 
such  that  the  corresponding  Lagrangian  relaxation(s) ,  if  solved,  would 
provide  "tighter"  bounds.  If  such  a  relaxation  does  not  possess  the 
integrality  property,  then  it  provides  an  equal  or  better  bound  compared 
to  that  from  an  LP  relaxation,  as  mentioned  in  Chapter  2, 

Two  alternative  formulations  of  problem  (P) ,  along  with  their 
Lagrangian  relaxations,  are  given  below. 


6.1.1  Alternative  Formulation  1 


Formulation  (API)  is  obtained  by  adding  the  constraints 
eik  Xij  -  Xk  *  for  all  i,j,  and  k  ,  to  problem  (P) ,  i.e,, 


(API) 


to 


X  £ 

i  j 

X 

i 

£  X 

i  j 


a  .  .  x  .  .  +  £  b,  y, 
ij  ij  k  k  k 

O) 

x.  .  =  1 
ij 

YJ 

(2) 

r  .  x  .  .  <  y. 

ljk  ij  -  k 

V’« 

(5') 

eik  xij  yk 

V  i,.i  ,k 

(52) 

xij  *  yk  =  0  or  1 

V  i ,  j  ,  k 

(6) 
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Since  *  1  or  0,  each  constraint  of  (52)  is  either  equivalent  to 

Xij  -  eik  ~  ^  °r  e^se  redundant  (if  *  0).  Problem  (API) 

thus  has,  at  ntost,  mnp  additional  constraints  relative  to  problem  (P)  . 

Two  Lagrangian  relaxations  are  now  considered  for  problem  (API). 

The  first  Lagrangian  relaxation  is  obtained  with  respect  to 
constraints  (5f)  by  introducing  nonnegative  Lagrange  multipliers 
v^  0  to  get 


Minimize 


£  Z  a  x..  +  Z  b  y  -  Z  v  (y  -Z  Z  r  ,  x.  \ 

i]  iJ  k  k  k  k  k  \*k  ±  .  ijk  ij ) 


subject  to 


(ALRlv) 


(2),  (52),  and  (6),  or 

equivalently, 

I  Z  x .  .  /  a  +  Z  v  r 

i  j  ^  l  lj  k  k 

ijk)  k  yk  (Vk  bk) 

(53) 

Z  x  .  =1 
i  ^ 

VJ 

(2) 

eikxij±  \ 

Vi, j ,k 

(52) 

xij  >yk  =  0orl 

Vi, j ,k 

(6) 

Another  Lagrangian  relaxation  of  problem  (API)  is  obtained  with 
respect  to  constraints  (51)  and  (52)  by  introducing  noni  egative  Lagrange 
multipliers  v  and  A  ,  respectively,  to  get 

K,  IJ  K. 


Minimize  Z  Z  a .  .  x .  .  -1-  Z  h  y. 

ij  1J  1J  k  k  k 


-  Z  v. 


(yk  - 1  *  r«k  xu) 


-  -  *tt  -u) 


subject  to  (2)  and  (6),  or  equivalently. 


Minimize 


(ALR1v,X)  I 


Subject  to 


V 


?,  E  x 

/a  +  E  v,  r .  .  _  +  Z  e  ,  X...  \ 

i  j  ij 

\  ij  k  k  ljk  k  :k  ljk/ 

-  E  \ 
k  k 

(54) 

Ex., 
i  1J 

=  1  VJ 

(2) 

xij  » 

yk  =  0  or  1  Vi,  j 

(6) 

For  this  problem,  the  solution  is: 


y 


k 


0 

1 


if 


if 


0  , 

0  , 


and  x .  . 
iJ 


1  if  i  minimizes 


l  ry* 


+  E 
k 


~£k 


over  £  . 


We  need  good  choices  of  Lagrange  multipliers  with  which  t;o 

solve  problem  (ALRl^) ,  and  of  Lagrange  multipliers  v^  and 

which  to  solve  problem  (ALR1  ,) .  Problem  (ALR1  )  d°es  not  possess  the 

V ,  A  V 

integrality  property,  thus  offering  the  hope  of  a  tigr.t  bound,  but  has 
more  constraints  and  is  difficult  to  solve  compared  to  problem  (ALR1 

v ,  A 

which,  on  the  other  hand,  involves  more  Lagrange  multipliers. 

6.1.2  Alternative  Formulation  2 

Another  formulation  of  problem  (P)  is  similar  to  problem  (API) 
except  for  a  modification  in  constraints  (51),  i.e.. 


i 

\  y 
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Minimize  Z  Z  a .  .  x  +  Z  b,y 
i  j  1J  lj  k  K  K 

Subject  to  £  x .  .  =  1 


(AP2) 


Z  Z  r .  ..  x.  .  <  1  Vk 

i  j  1Jk  «  ' 

eikxij  i  yk 


xij  ’  yk  =  0  or  1  Vi> j ,k 


A  Lagrangian  relaxation  with  respect  to  constraints  (55)  and  (52)  is 

obtained  by  introducing  nonnegative  Lagrange  multxpliers  v,  and  X  , 

K  1 J  K 

to  get 


Minimize  [  I  aj(  x..  +  Z  b  v 
,  .  ij  ij  .  kk 

i  j  J  J  k 


l'*  C - *«) 

v  (»k- eikxi3) 


Subject  to  (2)  and  (6),  or  equivalently, 


Minimize  I 


(ALR2  , ) ■ 
v,  A 


I  l  X.  . 

fa  .  +  £ 

v.  r  4-  Z  e 

i  j  1J 

l  to  k 

k  ijk  k 

+  Z  y, 

/  b,  -  Z  l 

xJ1(  \  -  r  v, 

k  K 

\k  fj 

ijkj  k  1 

E  x .  . 

=  1 

YJ 

i  1J 

V;  ’ 

yk  -  0  or 

1  Vi,  j 

For  this  problem,  the  solution  is: 


yk  =  0  if  EE  XiJt  ±  , 

‘  1  “  *  ‘ijk  i  \  • 
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and  x .  .  *  1  if  i  minimizes 
ij 


(a*j  +  l  \  r*Jk  +  l  e*k  XiJk)°Ver  f  ' 


Here  again,  we  need  good  choices  of  the  Lagrange  multipliers 

and  A..t  with  which  to  solve  problem  (ALR2  .). 
ijk  v  v ,  A 

6-1.3  Choice  of  Lagrange  Multipliers 

Each  of  the  relaxations  (ALRl  ,  )  and  (ALR2  ,)  involves  p  v. 

v,A  v,A  k 

Lagrange  multipliers  and  mnp  A...  multipliers.  If  we  have  good  choices 

i  j  k 

of  these  multipliers,  the  resulting  solutions  of  the  relaxed  problems 

should  provide  "tighter"  bounds  (because  of  the  additional  constraints) 

than  the  bound  from  relaxation  (LR  ).  Since  relaxaiions  (ALRl  «>  )  and 

v  v,A 

(ALR2  ,  )  are  similar  to  a  great  extent,  only  the  relaxation  (ALRl  ->) 

V  ,  A  V  >  A 

will  be  considered  for  further  discussion. 

By  looking  at  expression  (54)  of  the  formulation  (ALRl  x),  a 

v,  A 

meaningful  choice  of  the  Lagrange  multipliers  and  a  ^  appears 

to  follow  from  setting 


v.  +  Z  l  A,,,  «  b.  Vk 

k  i  .  ijk  k 

so  that  each  of  the  A...  can  be  chosen  as 

ijk 


bk  '  Vk 
^jk  n( Z  e.,  \ 


if  £ik  "  1  * 


0  otherwise 


The  solution  for  problem  (ALRl  .)  is  then  to  select,  from  each  column  j 

v,  A 


an  variable  which  minimizes  +  T  ^  +  E  ^  X^ 


over  £ 


Arbitrary  values  were  considered  for  the  (e.g.,  equal  to 

3/4  ,  1/2  b^  *  1/4  b^  *  and  0)  ,  the  A^^  were  then  computed  from 
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(58),  and  the  test  problem  with  three  designs  (m) ,  four  activities  (n) 
and  five  facilities  (p)  was  solved.  Three  cases  with  different  capacities 
(as  specified  in  Chapter  5,  Table  4)  were  tried  for  the  solutions  at 

the  initial  node.  The  results,  however,  were  not  conc'.usive  in  terms  of 
providing  a  meaningful  choice  of  the  Lagrange  multipliers  v^  (and  of 


Since  the  relaxation  (ALR1  x)  possesses  the  integrality  property, 

V,A 

a  choice  of  the  multipliers  as  the  optimal  values  of  the  dual  variables 
of  the  corresponding  linear  programming  problem  woulo  provide  a  solution 
as  good  as  the  LP  solution  (as  stated  in  Chapter  2) .  We  do  not  propose  to 
solve  linear  programs  as  a  part  of  our  branch-and -bound  methodology. 
However,  we  have  made  some  LP  runs,  basically  to  see  if  the  results  provide 
insight  leading  co  the  choice  of  the  Lagrange  multipliers,  and  also  to 
see  if  the  resulting  LP  solutions  are  "close”  to  tha  integer  solutions. 
These  results  are  given  below. 


The  LP  formulation  (API)  correspond ing  to  probxern  (API)  is: 

r 


Minimize 


Subject  to 


v 

The  constraints  x.. 

ij 


a 

ij 

x .  .  +  E  by 
ij  k  k  k 

(A) 

x .  . 

ij 

=  1 

VJ 

(2) 

rij  k 

x  .  .  <  y. 

ij  - 

V'c 

(5’) 

eik 

Xij  ±  yk 

Vij.k 

(52) 

M1 

Vk 

(15) 

x ,  . 

.  yk  >  0 

Vi , j  »k 

(16) 

<  1  are  implicit  in  constraints  (2). 


p=5 


Problem  (API)  was  solved  for  the  test  problem  with  m=3,  n=4 ,  and 

and  three  different  cases  for  the  capacities  (os  specified  in 
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Chapter  5,  Table  4)*  Each  case  was  solved  using  the  IMSL  (International 
Mathematical  and  Statistical  Library)  Code  ZX3LP  on  an  IBM  3031  at  The 
George  Washington  University. 

Note  that  the  formulation  (API)  has  up  to  mnp  more  constraints  than 
the  LP  formulation  (P)  given  in  Chapter  2.  For  our  test  problem,  this 
translates  into  solving  a  problem  of  17  variables  and  50  constraints 
corresponding  to  formulation  (API)  as  against  17  variables  and  14  con¬ 
straints  corresponding  to  formulation  (P) . 

Table  5  lists  the  solution  values  for  each  of  the  three  cases  with 
different  capacities  for  the  small  problem  with  thr^e  designs,  four  activi¬ 
ties  and  five  facilities.  The  solutions  to  problem*  (P)  and  (API), 
obtained  from  ZX3LP,  show  the  optimal  solution  values,  the  optimal  values 
of  the  variables  x_  and  y^  ,  and  the  optimal  values  of  the  dual 

variables  corresponding  to  the  Lagrangian  relaxations  (LR  )  and  (ALR1  >), 

V  V  ,  A 

i.e.,  v^  associated  with  the  capacity  constraints  (j!)  and 

associated  with  the  constraints  (52).  The  table  also  shows  Z(P),  and 
the  Lagrangian  solution  value  Z(LR^)  obtained  by  setting  v^  =  b^  for 

all  k  at  the  root  node,  i.e,,  Z(LR,  ) . 

b 


As  expected,  the  LP  solutions  for  each  of  the  three  cases  show 
Z(AP1)  to  be  considerably  higher  than  Z(P),  and  closer  t~  Z(P),  thereby 
providing  a  tighter  bound.  As  for  the  Lagrange  multipliers  v^  and 


X 


the  following  relationships  are  observed. 


z  z 

i  j 

A.  *  1, 
ijk 

1  bk 

V  k  ,  and 

vk  + 

Z  Z 

i  j 

Xijk  -  bk 

Vk  . 

Also , 

for 

Vk  = 

0  »  1  l.  Xijk  = 
l  J 

b  ,  and 
k 

for 

Vk± 

bK  •  \)  Xlik  ’ 

> 

o 

TABLE  5 

LP  AND  OTHER  SOLUTION  VALUES  FOR  A  TEST  PROBLEM 
(m=3,  n=4 ,  AND  p=5) 
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Although  the  relationship  among  various 


A. 

ijk 


values  is  not  apparent. 


the  above  observations  are  useful  in  further  exploring  some  good  choices 

of  the  Lagrange  multipliers  for  the  relaxation  (AI.R1  ,),  as  discussed 

v ,  A 


earlier . 


As  for  the  ’’closeness1’  of  the  LP  solution  to  that  of  the  integer 
solution,  most  of  the  solution  values  x  and  y  of  problem  (API) 

ij  k 

are  fractional,  and  their  rounding  off  to  0  or  1  does  not,  in  general, 
seem  to  correspond  to  the  optimal  integer  solution  values  .  and 

yk  of  problem  (P). 

Table  5  also  displays  Z(LR^)  at  the  root  node  for  each  of  the 
three  cases.  For  sk  =  3000  ,  Z(LRb>  =  Z(P)  ,  and  the  Lagrange  multi¬ 
pliers,  as  reflected  by  the  values  of  the  dual  variables  of  problem  (P) , 
are  equal  to  for  all  k  .  This  is  expected  from  Theorem  3  and  our 

discussion  of  the  integrality  property  in  Chapter  2.  Further,  the  dual 
variables  of  (P)  for  the  first  two  cases  (i.e.,  sR  -  7C0  Vk  ,  and 

s  =  400,...)  have  values  v  >  b.  from  Theorem  1. 


The  ZCLR^)  values  in  Table  5,  however,  take  no  consideration  of 
the  capacity  rule  and/or  the  facility  usage  rule  of  our  branch-and- 
bound  procedure.  These  rules,  by  fixing  certain  x_  values  to  0 

or  1,  and  by  forcing  certain  facilities  into  the  solution,  provide 

an  improved  lower  bound.  As  per  our  branch-and-bound  procedure  the 

improved  lower  bound  at  the  root  node  is  obtained  by  solving  problem 

(PR.).  For  example,  for  s.  =  700  Vk,  the  values  of  Z(LR,)  and  Z(PR, ) 
x  k.  b  1 

are  shown  in  Figure  8,  The  figure  also  shows  the  values  of  Z(P) , 

Z(APl),  and  Z(P).  The  branch-and-bound  procedure  rules  provide  an 

improved  value  of  the  lower  bound  Z(PR^)  compared  to  Z(P).  It  appears 

that  some  good  values  of  the  Lagrange  multipliers  of  the  relaxation 

(ALRl  .),  if  found,  could,  in  conjunction  with  these  rules,  provide 
V ,  A 

significant  improvement  over  Z(AP1),  and  without  the  need  to  solve  an 
LP  problem. 


j 
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6.2  Subgradient  Method 

It  was  mentioned  in  Chapter  2  that  setting  the  Lagrange  multipliers 
equal  to  for  all  k  provides  a  good  starting  point  in  solving 

the  Lagrangian  relaxation  (LRv)  of  problem  (P)  .  From  Theorem  3,  this 

choice  is  optima]  (in  terms  of  providing  the  tightest  lower  bound)  if 
the  resulting  solution  is  feasible  in  problem  (P) .  In  other  cases, 
i.e.,  where  the  resulting  solution  is  not  feasible  ir  problem  (P) ,  this 
choice  is  generally  not  optimal  and  it  is  possible  to  tighten  the  bounds 
by  considering  values  of  v^  b^  (from  Theorem  1) .  One  method  that 

seems  useful  in  providing  such  a  choice  is  the  subgradient  method. 

The  subgradient  method  is  an  adaptation  of  the  gradient  method 
in  which  gradients  are  replaced  by  subgradients.  Through  a  heuristic 
choice  of  the  step-size,  this  method  has  been  successfully  used  to 
provide  improved  bounds  and  sometimes  optimal  solutions  [for  details 
see  Held,  Wolfe,  and  Crowder  (1974),  Fisher  (1978),  and  Christofides 
(1980)].  The  fundamental  theoretical  result  is  that 

g 

Z(LR8)_»Z(D)  if  t8-*0  and  Z  tq-*°°  as  g  — ►  °°  , 

V  q=0 

where  t8  is  the  positive  step-size  t  for  the  gth.  iteration, 
and  Z( LR^8)  is  the  solution  value  of  the  relaxed  problem  (LR^)  using 

v^  values  obtained  at  the  gth  iteration. 

C’f'l 

In  the  case  of  problem  (P) ,  the  step-size  t&  for  iteration 
g  +  1,  given  that  we  have  a  solution  of  (LR^),  is  given  by 


X8+1[Z*  -  Z(LR  8)] 


"  J  j  rijk  Xij 


g|  |2 


(59) 


g+1 

where  is  a  scalar  and  generally  between  0  and  2,  and 

Z*  is  an  upper  bound  on  Z(LRv^),  frequently  obtained  by  applying 


a  heuristic  to  solve  problem  (P) . 
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Given  the  vector  of  multipliers 


,8+i 


is  generated  by 


II  r 

i  J 


ijk 


(60) 


where  we  enforce 


in  our  case  of  problem  (P)  (because  of 


Theorem  1) . 


Justification  for  these  rules  and  conputat ional  results  of 
applications  of  the  subgradient  method  are  given  in  Held  ett  ajL  (1974). 

The  scalar  A  is  generally  initiated  by  setting  A^  =  2  and  halving 
subsequent  A’s  whenever  the  resulting  solution  value  has  failed  to 
increase  in  some  fixed  number  of  iterations.  This  rule  has  performed 
well  empirically  [Held  et  al_  (1974)  and  Fisher  (1978)]. 

For  the  test  problem  with  three  designs,  four  activities,  five 
facilities,  and  the  capacities  =  400,  S2  *  400,  «*  1000, 

s^  =  400,  =  400,  the  Lagrangian  solution  obtained  it  the  root  node 

by  setting  =  b^_  for  all  k  ,  i,e,#  the  solution  to  problem  (PR^) 

is  infeasible  ir  problem  (P) ,  i.e.,  it  violates  the  capacity  constraints. 
It  seems  that  the  subgradient  method  could  be  useful  in  considering 
v  b  with  the  ultimate  purpose  of  obtaining  a  tighter  lower  bound, 

K  K 

and,  depending  on  the  revised  solution(s),  possible  improvement  in  the 
best  upper  bound.  Another  possiblity  could  be  to  first  arbitrarily 
increase  the  relevant  v^  values  by  a  small  percentage  of  the  b^ 

values  and  then  solve  problem  (LR^) >  hopefully  to  improve  the  lower 

bound;  and  thereafter  to  use  the  subgradient  method  for  obtaining 
subsequent  values  of  v^  ,  and  tightening  the  bounds. 


Both  ol  the  areas  discussed  above,  i.e.,  the  consideration  of 
alternative  f ormula t ions  of  problem  (P),  and  the  application  of  the 
subgradient  method,  and  their  combination,  seem  useful  for  continued 
research  in  terms  of  further  improving  the  branch-ani-bound  procedure 
for  solving  the.  multiactivity  multifacility  capacity-constrained  0-1 
assignment  problems. 
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APPENDIX  A 

ZIPCAP  LISTING  (REVISED) 


FORTRAN  IV  G  LEVEL  21 


MAIN 


DATE  =  80315 


11/07/24 


0001 


0002 

0003 

0004 

0005 


0006 

0007 


0000 

0009 

0010 

0011 

0012 

0013 


001<* 


0015 
0016 
OU  1  7 
0018 
0019 
0020 
U02  1 
0022 


C 

C 

C 

C 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

L 

c 

c 

c 

c 


c 

c 

c 


c 


c 


c 

c 

c 

c 


ZIPCAP,  A  ZERO-ONE  INTEGER  PROGRAM  IS  DESIGNED 
TO  SOLVE  MULTI  ACTIVITY  MUL T I F AC  I L I TY  CAPACITY- 
CONSTRAINED  PROBLEMS  HAVING  VARIABLE  AND  FIXED 
COSTS.  IT  ALSO  SOLVES  UNC At  AC  I TATED  PROBLEMS  AS  A 
SPECIAL  CASE 

INTEGER  0(35,35,30),  A<35,35),  CX(35,3f),  E ( 35 , 30 ) , 

1  0(30),  BSOLX ( 3  5  )  ,  0  SOL Y (30),  FL8 t 30 ) , F I X ( 35 ) , F I X I C 3 5  I  , 

2  FUBOO),  S  (  30  )  ,  SOLX(  35  )  ,  STX(122T> 

REAL  M I NC (35),  NMINC(35) 

DIMENSION  C( 35, 35 ) ,  DIFBR(35),  KT2C35),  MIND(3») 

INTEGER  BRO,  BR1 ,  FC,  FCU8 ,  P 

REAL  LOWB,  MAXDIF,  MTNSC 

********** ****QPTI ON S  AVAILABLE:  IINPT,  ICAPR,  I  STEP,  IUNCAP,£PS 
1 INP T  =  1  IF  INPUT  LISTING  DESIRED;  0  OTHERWISE 
I C APR  =1  IF  CAPACITY  RULE  TO  BE  USED;  0  OTHERWISE 
I STEP*0  IF  LISTING  OF  INTERMEDIATE  STEPS 
NOT  DESIRED.  IST£P=i  IF  SUMMARY  OF  BRANCH  L 
BOUND  NODES  DESIRED.  ISTEP=2  IF  DETAILED 
LISTING  OF  INTERMEDIATE  STE°S  DESIRED. 

I UNC AP  =  1  IF  SOLVING  AN  UNCAP AC  I TATED  PROBLEM, 

0  OTHERWISE. 

6PS=  A  FRACTIONAL  VALUE  IF  LUBOPT IMAL 
SOLUTION  DESIRED,  E.G.,  EPSILON  AS  0.005 
IMPLIES  SOLUTION  TO  BE  WITHIN  vC.5  PERCENT 
OF  THE  OPTIMAL  SOLUTION.  EPS=0«0  IF  OPTIMAL 
SOLUTION  DESIRED. 

E ELAPSED  TIME  IN  SECONDS,  IF  SPECIFIED,  AT 
WHICH  THE  NODE  AND  BOUNDS  RELATED  INFORMATION 
IS  PRINTED.  THIS  IS  USEFUL  IN  A  SITUATION  IF 
ISTEP=0  AND  THE  PROGRAM  TERMINATES  BEFORE 
REACHING  THE  FINAL  SOLUTION. 

INPUT  DATA********************************* 
READ  10,  IINPT,  ICAPR,  ISTEP,  IUNCAP,  EPS,  ET 
10  FORMAT  (411,  F6.5,  F10.3) 

M*  NUMBER  OF  DESIGNS 
N=  NUMBER  OF  ACTIVITIES 
P*  NUMBER  OF  FACILITIES 

READ  20,M,N,P 
20  FORMAT  (315) 

A (  I ,  J  )  :  VARIABLE  COST  MATRIX 
READ  30,  (IA|1,J),  I=1,M),J=1,N) 

30  FORMAT  (8110) 

BCKI :  FIXED  COST  V'  CTOR 
READ  30,  ( B ( K )  , K  =  l ,P) 

IF  (IUNCAP. EQ.l)  GO  TO  40 

S(K):  CAPACITY  LIMIT  VECTOR;  REQUIRED  ONLY 
IF  1 UNC A  P  =  0 
RE  AO  ,0,  ($(K),K*1,P> 

0(1, J,K):  CAPACITY  USAGE  MATRIX;  REQUIRED 
ONLY  IF  IUNCAP=0 

DO  32  K  =  1  , P 

READ  30  ,  ( ( DC  I  ,  J,K  )  ,  1*1 , M > , J= l ,N I 

32  CONTINUE 
DO  37  K  = 1 , P 
00  37  I  =  l  ,  M 

IF  (0(X,l,K).E0.O)  GO  TO  35 
t  (  1  ,K  )  s  I 
GO  TO  37 


00000010 

00000015 

00000020 

00000030 

00000040 

00000050 

00000060 

00000070 

00000080 

00000090 

ooooono 
00000120 
00000130 
00000140 
00000150 
00000160 
00000170 
00000180 
00000190 
00000200 
00000210 
00000220 
00000230 
00000240 
00000250 
00000253 
00000256 
00000260 
00000263 
00000266 
00000270 
00000273 
00000280 
00000290 
00000300 
00000310 
00000320 
00000330 
00000340 
00000350 
00000360 
00000370 
00000380 
00000390 
00000400 
00000410 
00000420 
0000 0430 
0000CH40 
00000450 
00000460 
00000470 
00000430 
000004^0 
OCOOCSOO 
000005 1 0 
00000520 
00000530 


AD-A102  583  GEORGE  WASHINGTON  UNIV  WASHINGTON  DC  PROGRAM  IN  LOGISTICS  F/G  12/1 

SOLVING  MULTIACTIVITY  MULTIFACILITY  CAPACITY-CONSTRAINED  0-1  AS— €TC(U> 
MAY  81  K  L  CHHABRA  N00014-B0-C-0169 

UNCLASSIFIED  SERIAL-T-441  NL 


END 
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0023 

35  E ( I  * K 1 *0 

00000540 

0024 

37  CONTINUE 

00000550 

0025 

GO  TO  90 

00000560 

C 

E 1 1 f K  ) :  DESIGN-FACILITY  MATRIX;  REQUIREO  ONLY 

00000570 

C 

IF  IUNCAPsI 

00000580 

0026 

*0  READ  45,ME(I,K),I*l,M>tK«l,P> 

00000590 

0027 

45  FORMAT  (8011) 

00000600 

0028 

DO  80  K*I , P 

00000610 

0029 

5  (  K  )  *N 

00000620 

0030 

00  75  I  *  l  ,  M 

00000630 

0031 

IF  (E(I,K).EQ.l>  GO  TO  65 

00000640 

0032 

00  60  JxlvN 

00000650 

0033 

Dili J,KI*0 

00000660 

0034 

60  CONTINUE 

00000670 

0035 

GO  TO  75 

00000680 

0036 

65  DO  70  Jxl,N 

00000690 

0037 

OC  I,J,K)*l 

00000700 

0038 

70  CONTINUE 

00000710 

0039 

75  CONTINUE 

00000720 

0040 

80  CONTINUE 

00000730 

c 

00000740 

0041 

90  PRINT  95*  IlNPTt  ICAPR*  ISTEP,  IUNCAP,  EPS,  ET 

00000750 

0042 

95  FORMAT  (  •  1*  t  •  OPTIONS  SELECTED  :  IINPT*',I1, 

00000760 

1  •  IC APR* # , 11 ,  •  ISTEP** ,11,  *  IUNCAP**, 11, 

00000770 

2  *  EPS** , F8 • 5,  *  ET*1,  F10.3///) 

00000780 

0043 

IF  (IINPT.EQ.O)  GO  TO  168 

00000790 

0044 

PRINT  100  »M , N , P 

00000800 

0045 

100  FORMAT  (*0*,  T55,  *  INPUT  DATA*, /IX,  T55,  f  — — —  — — *  ,//////lX  , 

00000810 

1T41,  'NUMBER  OF  DESIGNS  (M)**,  4X , I4//1X, T4l , 

00000820 

2 'NUMBER  OF  ACTIVITIES  IN  1  * • »  1X.I4//IX,  T<1» 

00000830 

3 • NUMBER  OF  FACILITIES  <P)**»IX,  14///) 

00000840 

0046 

PRINT  105 

00000850 

0047 

105  FORMAT  (  4X,  'VARIABLE  COST  MATRIX  A I  I , J ) • ,/4X, 

00000860 

1  * - »,/) 

00000870 

0048 

DO  110  1*1, M 

00000880 

0049 

110  PRINT  115,  I,  <A(I ,J),J*1,N> 

00000890 

00  50 

115  FORMAT  CO*,  T6,  •!**,  13,  4X,8I13,  4(/,  I4r,8113)) 

00000900 

0051 

PRINT  120 

00000910 

0052 

120  FORMAT ( *0* ,//4X, 'FIXED  COST  VECTOR  B(K)',/4X, 

00000920 

1* - - '  ,  /) 

00000930 

0053 

PRINT  122,  (B(K) ,K*1 ,P ) 

00000940 

0054 

1 22  FORMAT  ('O',  T15,  8113,  31/,  I4X, 8I13II 

00000950 

0055 

PRINT  125 

00000960 

0056 

125  F0RMAT('0',//4X, 'CAPACITY  LIMIT  VECTOR  S(M*,/4X, 

00000970 

1« - -  ■  -»  ,  /| 

00000980 

0057 

PRINT  128 1  (  S !  K ) ,  K*  1  *  P 1 

00000990 

0058 

128  FORMAT  I'O',  T15,  3113,  31/,  14X, 8113)1 

00001000 

0059 

PRINT  130 

00001010 

0060 

130  FORMAT! 'O' ,//4X,  'CAPACITY  USAGE  MATRIX  0 ( I , J ,K • • , /4X, 

00001020 

1'-  - - ',/) 

00001030 

0061 

DO  150  K*i,P 

00001040 

0062 

PRINT  135, K 

00001050 

0063 

135  FORMAT  ( *0  * *//5X , •*  =  • , 13/ 1 

00001060 

Q0o4 

00  145  1=1, M 

OCOO 1070 

0065 

PRINT  140,1  ,(D( I,J,K) ,  J*1 ,N  ) 

00001080 

0066 

l 40  FORMAT  ('O',  T6,  '1=',  13,  4X,8I13,  4(/,  l4X,8113>) 

00001090 

00o7 

145  CONTINUE 

00001 100 

0068 

150  CONTINUE 

GGOOlilO 
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t 


0069 

PRINT  155 

00001120 

0070 

155 

FORMATt *0* ,//4X, •DESIGN-FACILITY  MATRIX  E(I»K!*«/4X, 

00001130 

( 

!•■■■  ■■■■  •■■■■■  - — ',/) 

00001140 

0071 

00  160  I* 1 »  M 

00001150 

0072 

PRINT  158,  I,  (Eli  ,tOtK*ltP> 

00001160 

( 

0073 

150 

FORMAT  ('O',  T6,  •!=',  13,  4X,8U3,  3(/,  14X,8113>) 

00001170 

0074 

loO 

CON . INUE 

00001180 

0075 

i6e 

IF  /ISTEP.EO.O)  GO  TO  190 

00001190 

3 

0076 

IF  (ISTEP.EQ.il  GO  TO  175 

00001200 

0077 

PRINT  170 

00001210 

0078 

170 

FORMAT  < 'O', ///55X, 'DETAILED  LISTING  OF  HEPS',/) 

00001220 

,  \ 

0079 

GO  TO  190 

00001230 

0080 

175 

PRINT  180 

00001240 

0081 

180 

FORMAT  ( 'O', ///55X, 'SUMMARY  OF  STEPS',/) 

00001250 

o 

C 

00001260 

0082 

190 

BUB*9999999. 

00001270 

0083 

BUBS  =  BUB/  ( 1.0+EPS ) 

00001280 

o 

0084 

NSX*0 

00001290 

0085 

NOO=1 

00001310 

0086 

I BNOD  =  1 

00001315 

c 

0087 

INE  T  =  0 

00001320 

0088 

I NS  ET  =0 

00001330 

0089 

DO  205  J  =  l  ,N 

00001390 

c 

0090 

f: X( J)=0 

00001400 

0091 

KT2< J)=0 

00001410 

0092 

DO  205  1=1, M 

00001420 

c 

0093 

CX(1,J)=0 

00001430 

0094 

205 

CONTINUE 

00001433 

0095 

LQ  1  =0 

00001436 

(j. 

0096 

LG  2=0 

00001440 

0097 

LR2=0 

00001443 

0098 

CALL  TIMET ( ITO) 

000014*5 

0 

0099 

IF  (ISTEP.EQ.O)  GO  TO  208 

00001448 

0100 

PRINT  220, NOD 

00001450 

0101 

208 

IF (NSX.EQ.OI  GO  TO  283 

00001453 

c 

c 

CXU,J)  CONTAINS  FIXED  ANO  FREE  X(I,J>  VARIABLES. 

00001456 

c 

STX(INS)  CONTAINS  FIXED  X(I,J)  VARIABLES. 

00001460 

c 

CX(I,J)  AND  STX(INS)  ARE  UPDATFD  BY  THE  CAPACITY 

00001480 

o 

c 

RULE,  THE  BOUNDING  RULE,  AND  THE  RULE  FOR 

00001490 

c 

BRANCHING  AND  BACKTRACKING. 

00001500 

c 

IN  CX(I,JI  A  FIXED  VARIABLE  ,  S  RECORDED  AS  1  OR 

00001505 

o 

c 

2,  ANO  A  FREE  VARIABLE  AS  0. 

00001510 

c 

A  VALUE  OF  l  IMPLIES  THAT  THAT  PARTICULAR  VARIABLE 

00001515 

c 

IS  FIXED,  AND  FIX(J)  IS  SET  EQUAL  TO  I  IMPLYING 

00001520 

o 

c 

THAT  COLUMN  J  HAS  A  FIXED  VARIABLE  OF  VALUE  I. 

00001525 

c 

A  VALUE  OF  2  IMPLIES  THAT  THAT  PARTICULAR  VARIABLE 

00001530 

c 

SHOULD  NOT  BE  CONSIDERED  FOR  CURRENT  COMPUTATIONS. 

00001535 

0 

c 

AN  X(I,J)  RECORDED  IN  CXU,J)  A$  1  DUE  TO  THE 

000015*0 

c 

BRANCHING  RULE  IS  RECORDtD  IN  STXf.’NS)  AS  X*100  +  J. 

00001545 

c 

AN  X(I,J>  RECORDED  IN  CX(I,J>  A$  l  DUE  TO  THE 

00001550 

i 

c 

CAPACITY  RULE  OR  THE  BOUNDING  Rl,LE  IS  RECORDED  IN 

00001555 

c 

S  T  X  (  INS  )  AS  (XMOO*J)MOOCOOO. 

00001560 

c 

AN  RECORDED  IN  CX(I,J)  AS  2  IS  RECORDED  IN 

0000 1565 

c 

ST  X ( I  NS )  AS  -( X*l GO+J  1-1000000 . 

00001570 

0102 

210 

IF  (l$TEP.tQ.O>  GO  TO  225 

00001580 

0103 

215 

PRINT  220, NOD 

00001590 

0104 

220 

FORMAT  (  •0,,//6X,*N0DE  NUMBER',  17/1 

00001600 

c 

00001610 
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C 

o 

© 

o 

o 

o 

o 

c 

c 

o 

o 

o 


C 

BRO  IS  THE  RIGHT  BRANCHING  VARIABLE 

00001615 

0105 

225 

LX«BR0 

00001620 

0106 

I X*LX/ 100 

00001630 

0107 

JX*LX~IX*10Q 

00001640 

0108 

CX( IX»JX)*2 

00001650 

0109 

KT2<  JX)*KT2<  JXW1 

00001660 

0110 

F IX ( JX ) *0 

00001720 

0111 

L01*LQl-l 

00001725 

0112 

IF  UT2I  JX».LT.CM-1>>  GO  TO  270 

00001730 

0113 

00  255  1*1 

00001740 

0114 

IF  (CX«  I,JX>.EQ.2>  GO  TO  255 

00001750 

0115 

CXII  ,JXI*i 

00001760 

0116 

NSX*NSX*1 

00001763 

0117 

STX ( NSX 1*  (  I*1004JX>41000000 

00001766 

0118 

FI  X ( JX ) *1 

00001770 

0119 

LQ 1*LQ1>1 

00001780 

0120 

FI  XI ( JX 1*1 

00001790 

0121 

GO  TO  270 

0C001800 

0122 

255 

CONTINUE 

00001810 

0123 

270 

LQ2*0 

00001820 

0124 

LR2*0 

00001825 

0125 

GO  TO  233 

00001830 

0126 

272 

IF  ( I  STEP* EC . 0 )  GO  TO  276 

00001840 

0127 

PRINT  220.NOO 

00001850 

c 

00001853 

c 

BR1  IS  THE  LEFT  BRANCHING  VARIABLE 

00001856 

0128 

276 

LQ2*0 

00001860 

0129 

LR2*0 

00001866 

0130 

LX*BR1 

00001870 

0131 

1X*LX/100 

00001875 

0132 

JX*LX-IX*lOO 

00001860 

0133 

CX<IX»JX)*1 

00001885 

0134 

FIX ( JX )*1 

00001890 

0135 

LQ1*LQU1 

00001892 

0136 

00  279  I  =  l*M 

00001895 

0137 

IF  *IX.EQ.I>  GO  TO  281 

00001897 

0138 

279 

CONTINUE 

00001900 

0139 

281 

FIXII  JXl*lX 

00001902 

0140 

283 

IF  ( 1  STEP .NE «2 )  GO  TO  303 

00001905 

0141 

285 

DO  295  1*1 ,M 

00001910 

0142 

PRINT  290  y  I  « (CX( I v J ) t J*1 «  N) 

00001920 

0143 

290 

FORMAT  l/5X,'CX(  I,Jt*  ,*X,’1>*»13« 

2  X  v  20 14/23X  y  20141 

00001930 

0144 

295 

CONTINUE 

00001940 

0145 

PRINT  297,  < F I X 1 J), J»l,Nt 

00001950 

0146 

297 

FORMAT  </5X, *FIX ( J)»y 12X»  20I4/23X, 

20141 

00001960 

c 

00001970 

c 

AND  UPDATE  CXI  I  ,J)  AND 

STXC INSt . 

00001980 

0147 

303 

00  307  K* l yP 

00002000 

0148 

FLB ( K ) =0 

00002015 

0149 

307 

CONTINUE 

00002025 

0150 

310 

00  2000  K*lyP 

00002030 

c 

F1N0  THE  SUM  OF  MINIMUM  OllyJ.X)  OVER 

EACH  J  FOR  A 

000U2040 

c 

GIVEN  Ky  I.E.*  MINSC* 

SUM  OF  v I  NO  1  J  ) 

00002050 

0151 

M I N$0*0 

000020O0 

0152 

00  500  J*1 »N 

00002070 

0153 

IFIFIXl  Jl.EQ.O)  GO  TO  350 

00002030 

c 

IF  FlX(J)*ly  SET  MINCHJMOlIyJyKl  FOR 

CXI  I «  J 1  *  1 

00002090 

c 

ANO  MOVE  TO  NEXT  COLUMN  J 

00002 100 

* 


r 
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0154 

I NDI=F I X I  (  J  ) 

00002110 

0155 

HINDI J>=0(  IN01 ,  J,M 

00002120 

< 

0156 

GO 

TO  800 

00002130 

0157 

350 

LK  *0 

00002160 

0158 

1*1 

00002170 

(. 

0159 

HINDI J)*D(I, J,K) 

00002180 

C 

SKIP  DUtJtK)  WHEN  CX ( I ,J )*2  C  HOVE  TO  NEXT  ROW  I 

00002190 

OloO 

400 

IF (CX ( I ,J) .EQ.2 1  GO  TO  600 

00002200 

0161 

500 

TF(0(I, J,K1.LT.MIN0(  Jll  HINDU  )*0(  It  JrX  ) 

00002210 

01o2 

GO  TO  700 

00002220 

0163 

600 

LK*LK*1 

00002230 

0 

0164 

IFU.GT.LKI  GO  TO  700 

00002240 

0165 

i*m 

00002250 

0166 

HINDI Jl=0< I,J,K1 

00002260 

O 

0167 

GO  TO  750 

00002270 

0168 

700 

1*1*1 

00002280 

0169 

750 

IF ( I • LE  «H 1  GO  TO  400 

00002290 

C 

0170 

800 

M1NSD=MINSD*MIND(J> 

00002300 

0171 

900 

CONTINUE 

00002310 

0172 

910 

IF 

1 I STEP.NE .2)  GO  TO  960 

00002320 

o 

0173 

PRINT  950f  Kv  MlNSOy (M1ND(J) ,  J  *1 ,N ) 

00002330 

0174 

950 

FORMAT  ( • 0  * , #  K, M INSO , (MIND ( J ) » J*1 ,N 1 • , 10I10»4( /f44X, 81 101 ) 

00002340 

0175 

960 

IF  ( MI NSO • EG .0 )  GO  TO  975 

00002342 

o 

0176 

965 

IF  (  FLBUI.EQ.l)  GO  TO  975 

00002344 

0177 

970 

F  .BCKl*l 

00002346 

0178 

975 

IF  (IUNCAP.EO.l)  GO  TO  2000 

00002348 

o 

0179 

978 

IF  (ICAPR.EQ.O)  GO  TO  2000 

00002349 

c 

FIND  BALANCE  AVAILABLE  CAPACITY  I8AL0  FOR  A  GIVEN  K 

00002350 

c 

IF  IBALD  IS  NEGATIVE,  THEN  BACKTRACK. 

00002360 

L 

0130 

980 

IBALD*S(K1-MINS0 

00002380 

0181 

1000 

IF 

(IBALD. LT. 01  GO  TO  6200 

00002390 

0182 

00 

1500  J*i,N 

00002400 

0 

c 

SKIP  COLUMN  J  IF  FlXU!*i 

00002410 

0183 

IF  CFIXUI.EQ.U  GO  TO  1500 

00002420 

0134 

DO 

1300  1*1, M 

00002430 

c 

c 

SKIP  ROW  I  IF  CX ( I , J 1 *2 

00002440 

0185 

1100 

IF(CX( I,J1.E0.2I  GO  TO  1300 

00002450 

c 

COMPUTE  DIFFERENCE  BETWEEN  D(I«J,K)  AND  MINO(J). 

00002470 

o 

c 

IF  IT  IS  MORE  THAN  AVA1LABE  BALANCE,  SET  CX(I,JI*2 

00002480 

0186 

1200 

1DIFG=D(I  ,  4*  K  )«-M  IND  ( J  ) 

00002490 

0197 

IF  (IIDIFD-JBALD).LE.OI  GO  TO  1300 

00002510 

o 

0188 

CX  r 1 , J ) =2 

00002520 

0139 

NS  X=NSX  *1 

00002523 

0190 

STX ( NSX )=— ( I* lOO* J )— 1000000 

00002526 

c 

LQ2  COUNTS  THE  NUMBER  OF  CX(I,J)  VALUES  SET  EQUAL 

00002530 

c 

TO  2  IN  A  CYCLE 

00002540 

0191 

LQ2=LQ2  *1 

00002550 

o 

c 

KT2 ( J )  KEEPS  AN  ACCOUNT  OF  CX ( I , J 1  VALUES  SET  EQUAL 

00002560 

c 

TO  2  FOR  COLUMN  J 

00002570 

0192 

KT2 ( J ) =KT2 ( J 1 *1 

00002580 

< 

c 

FOR  COLUMN  J,  IF  ALL  BUT  ONE  CX(I,J)  VALUES  ARE 

00002590 

c 

EQUAL  TO  2,  SET  THAT  CXf!,Jl*l  C  SET  F1XCJ1M 

00002600 

0193 

IF (KT2< J) .LT. (M-ll )  GO  TO  1300 

00002610 

1 

0194 

00 

1250  LR=l,M 

00002620 

0195 

IF(CXCLR,J».EQ.2>  GO  TO  1250 

00002630 

0196 

C a  ('  R,J  1  =  1 

00002640 

0197 

NS  X=NS  X ♦ l 

00002643 

0198 

STX (NSX 1=  <LR* 100*31*1000000 

00002646 

89 
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0199 

FI  X ( J  )  *  l 

00002650 

c 

LOl  KEEPS  AN  ACCOUNT  OF  COLUMNS  FOR  WHICH  FIX(J)*t 

00002655 

0200 

LQl«LQi+l 

00002660 

c 

FIXIIJI  SPECIFIES  INDEX  I  FOR  WHICH  FIX<J)«I 

00002662 

0201 

FIXK  J)*LR 

00002665 

0202 

GO  TO  1500 

00002670 

0203 

1250 

CONTINUE 

00007680 

0204 

1300 

CONTINUE 

00002690 

0205 

1500 

CONTINUE 

00002700 

0206 

1800 

IF  ( 1STEP.NE.2)  60  TO  2000 

00002710 

0207 

PRINT  1900 •  K*  LQ2 1  LQ1 

00002720 

0208 

1900 

FORMAT  ( *0  •  »  *K* • *  13*  *  L02='»I3,  •  LQ1='»  13) 

00002730 

0209 

DO  1930  I*l*M 

00002740 

0210 

PRINT  290*  I*(CX(I*J)*J=I*N> 

00002750 

0211 

1930 

CONTINUE 

00002770 

0212 

PRINT  297,  (FIX(JI»J*l»NI 

00002780 

0213 

2000 

CONTINUE 

00002800 

C 

A  CYCLE  EXAMINES  ALL  THE  FACILITIES. 

00002803 

C 

IF  IN  A  CYCLE,  THE  CAPACITY  RULE  RESULTS  IN  SETTING 

00002810 

C 

ADDITIONAL  CX(I,J)  VALUES  EQUAL  TO  2,  THEN  REPEAT 

00002820 

c 

THE  CYCLE.  BUT  IF  FIX|J)«1  FOR  ALt  J,  THEN  DO  NOT 

00002830 

c 

REPEAT  THE  CYCLE. 

00002835 

0214 

IF  (LQl.EQ.N)  GO  TO  2400 

00002840 

0215 

IF  (LQ2.EQ.LR2)  GO  TO  2400 

00002845 

0216 

2200 

LR2SLQ2 

00002860 

0217 

GO  TO  310 

00002870 

C 

00002880 

C 

UPDATE  VECTOR  OF  FACILITIES  FL3CK)  FOR  COMPUTING 

00002890 

c 

C(ItJ)  MATRIX  C  LOWER  SOUND.  IT  HAS  VALUE  l  IF  A 

00002900 

c 

FACILITY  IS  USEDt  OTHERWISE  IT  HAS  0  VALUE. 

00002910 

G218 

2400 

DO  3000  J*l,N 

00002950 

0219 

IF  (FIX(J).EQ.O)  GO  TO  3000 

00002960 

0220 

INDI *F I XI ( J ) 

00002970 

0221 

DO  2550  K *  1  * P 

00002990 

0222 

IF  (EUNDltKKEQ.O)  GO  TO  2550 

00003000 

0223 

IF  ( FL8 ( K ) • EQ. 1 )  GO  TO  2550 

00003010 

0224 

FLB(K)=l 

00003020 

0225 

2550 

CONTINUE 

000030  30 

0226 

3000 

CONTINUE 

00003060 

0227 

IF  USTEP.N6.2)  GO  TO  3150 

00003070 

0228 

PRINT  3100,  (FLB(K)  ,K=l,P) 

00003080 

0229 

3100 

FORMA  T<  •0#»MFLB  IK)  ,K*ltPI  *  ,  2014/lbX, 2014 ) 

00003090 

C 

COMPUTE  COST  MATRIX  CU.J)  FOR  THE  RELAXED  PROBLEM 

00003100 

0230 

3150 

DO  3400  J*I,N 

00003110 

0231 

DO  3300  I *  I , M 

00003120 

02  32 

BSUMeO.O 

00003130 

0233 

DO  3200  K  * 1 , P 

00003140 

0234 

IF  (FLB(K) .60.1 )  GO  TO  3200 

00003150 

0235 

IF  <EU,K).EQ.O)  GO  TO  3200 

00003160 

0236 

B$UM«8SUM+<a(K)  *  IFLOATtOf I* JtK))/  FLOAT(SUM)) 

00003170 

0237 

3200 

CONTINUE 

00003180 

0238 

3250 

Cl  If J)*A( It J)4BSUM 

00003190 

02  39 

3300 

CONTINUE 

00003200 

0240 

3400 

CONTINUE 

00003210 

0241 

IF  IISTEP.NE.2)  GO  TO  3445 

00003220 

0242 

DO  3430  1*1, M 

00003230 

0243 

PRINT  3420,  I,  (C(I«J),J=l,N) 

00003250 

0244 

3420 

FORMAT  (/5X ,  •C(I,J)*,5X,  •!«•,  13, 2X,  5F15.4, 

00003260 
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1  6</23X,  5F15.4II 

00003265 

0245 

3430 

CONTINUE 

00003270 

C 

FIND  SUM  OF  MINIMUM  C<I,J>  VALUES  OVER 

EACH  J* 

00003290 

C 

I  *£•  f  MINSOSUM  OF  MINCU). 

00003300 

C 

IF  F  I  X(  J  )  *1 »  THEN  MINC(  jl^CU,  J)  WHERE 

CXII,J)«1 

00003310 

0  246 

3445 

MINSC=0.0 

00003320 

0247 

DO  3900  J  *  1 » N 

00003340 

0240 

IF  (FIX(J).EQ.O)  GO  TO  3500 

00003350 

0249 

INDIsFIXK  J) 

00003360 

0250 

MINC ( J 1 -C ( INDI«J) 

00003370 

0251 

S0LX( J J*INDI 

00003380 

0252 

GO  TO  3850 

00003410 

0253 

3500 

LK=0 

00003430 

0254 

i  =  l 

00003440 

C 

SKIP  C(I,J)  ELEMENT  IF  CX(y,J)=2  t  MOVE 

TO  NEXT  I 

00003470 

0255 

3550 

IF  (CXU»J).EQ.2)  GO  TO  370: 

00003480 

0256 

IF  ((l-LK).EO.l)  GO  TO  3600 

00003A85 

0257 

IF  (C(ItJ).GE.MlNC(J) )  GO  TO  3750 

00003490 

0258 

3600 

MlNC(J)*CI I,J> 

00003500 

02  59 

IMINrl 

00003510 

0260 

GO  TO  3750 

00003520 

0261 

3700 

LK=LK*l 

00003530 

02  62 

3750 

I*Ul 

00003590 

0263 

3800 

IF  (I »L  £*  M )  GO  TO  3550 

00003600 

02o4 

SOLX(J)*IMIN 

00003610 

0265 

3850 

MINSC=M 1NSC*M INC  1 J  J 

00003620 

0266 

3900 

CONTINUE 

00003630 

026? 

IF  USTEP.NE.2I  GO  TO  3940 

00003640 

02o8 

DO  3720  J= 1 » N 

00003650 

0269 

PRINT  3910  »  J,M1NC( J) ,S0LX( J) 

00003660 

0270 

3910 

FORMAT  (•OS*J*MINCU)*SOLXtJ>#,  I5tF15.4,16) 

00003670 

0271 

3920 

CONTINUE 

00003680 

C 

COMPUTE  FIXED  COST  FC  FOR  L( WER  BOUND 

00003710 

0272 

3940 

FC=0 

00003720 

0273 

DO  4000  Ka 1 f P 

00003730 

02  74 

IF  {FLB(K).EQ.O)  GO  TO  4000 

00003740 

02  75 

3950 

FC=FC*B(K) 

00003750 

0276 

'  000 

CONTINUE 

00003760 

c 

00003770 

0277 

4050 

LOWB=MINSC»FC 

00003780 

0278 

IF  (ISTEP.fcQ.OJ  GO  TO  4150 

00003790 

0279 

PRINT  4120,  MINSC,  FC ,  LOWB 

00003800 

0280 

4120 

FORMAT  (*0't*  MINSC,  FC,  LOWB  *,  F15.4,  IIS,  F15.4J 

00003610 

C 

COMPARE  LOWER  BOUND  WITH  BEST  UPPER  BOUND  STAR 

00003820 

C 

BUBS  WHICH  EQUALS  BU8/ll*EPS>.  IF  LOWB 

IS 

00003830 

c 

GREATER  THAN  OR  EQUAL  TO  BUBSt  THEN  BACKTRACK 

00003840 

0281 

4150 

IF  ( LOWB .GE. SUBS )G0  TO  6200 

00003850 

c 

CHECK  IF  CURRENT  SOLUTION  SATISFIES  CAPACITY 

00003880 

c 

constraints 

00003890 

0202 

4200 

IF  CIUNCAP.EC.1)  GO  TO  4420 

00003900 

0203 

4210 

DO  4400  K= 1  * P 

00003910 

0284 

NSUMOsO 

00003920 

0285 

00  4300  J*  1  *N 

00003930 

C20o 

1X*SCLX ( J ) 

00003950 

0287 

N$UMD*NSUMD>DtIXf J,  K) 

00003960 

0268 

o 

o 

■* 

CONTINUE 

00003970 

0289 

IF  (IST6P.Nfc.2J  GO  TO  4320 

00003900 

0290 

PRlNt  4310*  KfNSUMD 

00003990 
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0291 

0292 

0293 

0294 


0295 

0296 

0297 

0290 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 

C312 

0313 

0314 

0315 

0316 


0317 

0318 

0319 

0320 

0321 

0322 

0323 

0324 

0325 

032o 


0327 

0328 


4310  FORMAT  («0»,  9K,NSUM0 • #2I 101 

4320  IF(N$UM0.LE.S<K) )  GO  TO  4400 

GO  TO  5100 
4400  CONTINUE 

C  ♦•••♦••♦•♦♦♦•♦COMPUTE  UPPER  BOUND  UPB  IF  CAPACITY  CONSTRAINTS 

C  ARE  SATISFIED. 

C  UPBeSUM  OF  A( I , J ) +F IXED  COST  FCU3  BASED  ON 

C  SOLUTION  VECTOR  SOLX(J) 

C  VECTOR  OF  FACILITIES  FOR  UPPER  BOUND  FUB(K)  MAS 

C  VALUES  1  OR  0  BASED  ON  FACI1ILY  USED  OR  OTHERWISE 

4420  DO  4450  K»1,P 
FU8(K 1*0 
4450  CONTINUE 


NSUMA^O 

FCUB*0 

4500  DO  4650  J*1,N 
IX  =SOLX (J) 

NSUMA*NSUMA4A|IXt J) 

4550  00  4600  K*l,P 

IF ( E ( IX ,K ) «E0 .0  )  GO  TO  4600 
IF(FUB(K).EQ.l)  GO  TO  4600 
FUB(K)»1 
FCUB* FCUB ♦B ( K I 
4600  CONTINUE 
4650  CONTINUE 

IF  (ISTEP.NE.2J  GO  TO  4700 
PRINT  4o60,  (FUB(K) ,K=1, P) 

4660  FQRMAT( »0# *• (FUB(K) ,K=1,P)  •,  2014/ 16X , 20 14) 

4700  UPB*NSUMA*FCUB 

4708  IF  (1STEP.EQ.0)  GO  TO  4750 

PRINT  47 10 v  NSUMA,  FCUB,  UPBt  BUB,  BUBS 
4710  FORMAT ( *0# , • NSUMA ,  FCUB,  UPB,  BUB,  BUBS  •,2110,  2F15.4I 
C  COMPARE  UPPER  BOUND  WITH  BEST  UPPER  BOUND 

C  IF  UPB  IS  LESS  THAN  BUB,  SET  IT  AS  BUB  ANO 

C  NOTE  THE  SOLUTION 

4750  IF  (UPB. GE. BUB)  GO  TO  5100 
4770  BUB*UPB 

BUBS*  BUB/  (l.O^EPS) 

IBNOO=NOO 

PRINT  4780,  IBNOD,  BUB,  BUBS 
4780  FORMAT  (,0‘,  • IBNOD ,  BUB,  BUBS9 ,  110,  2F15.4) 

DO  4800  J* l ,N 
4800  BSOLXC J)*SOLX( J) 

00  4850  K*  1 , P 
4050  BSOLY(K)»FUB(K) 


C  ••••••••••••••COMPARE  LOWB  WITH  BUBS.  IF  LOWB  IS  GREATER 

C  THAN  OR  EQUAL  TO  BUBS,  THEN  BACK  TRACK 

4900  IF  (LOWB.GE.BUBS)GO  TO  6200 

C  ♦♦♦♦♦♦•♦♦•♦♦♦♦  IF  FIXU)  VALUES  ARE  1  FOR  EACH  J,  THEN  BACKTRACK 

5100  IF  (LQI.EO.N)  GO  TO  6200 

C  »»»**+*+»+ ,+»*AppLy  THE  BOUNDING  RUL £•♦•••••♦••*• * ••*♦•••♦♦ •♦* 

C  IF  THE  OlFFLRENCE  BETWEEN  C(I,J>  AND  MINC(J)  IS 

C  GREATER  THAN  THE  DIFFERENCE  BE TUtfcN  BUBS  AND 

C  LOWB,  THEN  CX (  I  ,  J  )  =  2 

C  •♦••••••♦•••♦♦APPLY  BRANCHING  RULE  ANO  FIND  BRl,  THE  NFXT 

C  VARIABLE  FOR  LEFT  BRANCHING. 

C  FINO  NMINC(J),  THE  NEXT  HIGHER  VALUE  THAN  MINC(J) 

C  ANO  OIFBR(J),  THE  DIFFERENCE  BETWEEN  THEM. 


00004000 
00004010 
00004020 
00004030 
00004040 
00004050 
00004060 
00004070 
00004080 
00004090 
00004100 
00004110 
00004120 
00004130 
00004140 
OOC  4150 
00004170 
00004190 
00004190 
00004200 
00004210 
00004220 
00004230 
00004240 
00004250 
00004260 
00004270 
00004280 
00004290 
00004300 
00004310 
00004320 
00004330 
00004340 
00004350 
00004360 
00004370 
00004380 
00004385 
00004386 
00004388 
00004390 
00004400 
00004410 
00004420 
00004430 
00004440 
00004450 
00004480 
00004500 
00004510 
00004515 
00004520 
00004525 
00004530 
00004540 
00004550 
0000*555 
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0383 

5810 

FORMAT  < 1 0 • »  * Jt NMINC ( J  )  » MINC (J)*DIFBR(J) #  ,  I5.3F15.4) 

00004963 

0384 

5820 

CONTINUE 

00004966 

C 

00004970 

0385 

5840 

DO  5900  J*1,N 

00004980 

0386 

IF  (J.NE.LJ)  GO  TO  5900 

00004990 

0387 

5050 

BRl=SOLX( J)*100*J 

00005000 

0388 

IF  (1STEP.EQ.0)  GO  TO  6020 

00005010 

0339 

PRINT  5880,  8R1 

00005020 

0390 

5880 

FORMAT (  'OS*  BR1M10) 

00005030 

0391 

GO  TO  6020 

00005040 

0392 

5900 

CONTINUE 

00005050 

C 

00005060 

c 

NSX  REPRESENTS  THE  NUMBER  OF  VARIABLES  IN  STXIINS) 

00005070 

0393 

6020 

NSX=NSXU 

00005090 

0394 

6040 

STX ( NSX ) *BR1 

00005100 

0395 

IF  (ISTEP.NE.2)  GO  TO  6100 

00005150 

0396 

PRINT  6088,  (  STX  ( INS  1  *  INS  =  l,NSXl 

00005160 

0397 

o088 

FORMA  T  (  •  0  •  ,  •  STXIINS)1,  10110,  122</,  lOX.lOIlOM 

00005170 

C 

**«*****»*****M0vE  TO  THE  NEXT  (LEFT  BRANCH)  NODS  AND  APPLY 

00005220 

C 

CAPACITY  RULE 

00005230 

0398 

6100 

NOD=NOD*l 

00005240 

0399 

6110 

IF  (ET.EQ.O.O)  GO  TO  6150 

0000524 2 

0400 

IFUNSET.EO.il  GO  TO  6147 

00005244 

0401 

IF  (INET.EG.l)  GO  TO  6150 

00005246 

0402 

CALL  TIMETCINTI 

00005248 

0403 

ELTN=( INT-II0)*26.04E-6 

00005250 

0404 

IF  (ELTN.LT.ET)  GO  TO  6150 

00005253 

0405 

6120 

PRINT  6125,  NOD,  ELTN,  BUB,  BUBS,  IBNOD 

00005256 

04  06 

6125 

FORMAT  COS  'WAS  AT  NODE1, 16,  •  AT  ELAPSED  TIME  *•,  FI0.4, 

00005260 

1 

•  SECQNOS. • ,/lX,  •  BUB=#,F15.4,  •  BUB $* * , F I 5 . 4 , 

00005263 

2 

•  AT  NODE= • , I 7 ) 

00005266 

0407 

I B  UB  XBUB 

00005267 

0408 

IF  (IBUB.EQ. 9999999)  GO  TO  6146 

00005268 

0409 

6130 

PRINT  6135,  (BSOLXi J),J=1,N) 

00005270 

0410 

6135 

FORMAT ( *0 • ,  "SOLUTION  CORRESPONDING  TO  BUB  IS*,  //IX, 

00005273 

1 

* (BSOLX< J),  J  =  i,N)  •, 10 18,3 (/18X, 1018) ) 

00005276 

0411 

6140 

PRINT  6145,  (BSOLY(K),  K=I,P) 

00005280 

0412 

6145 

F0RMAT(/1X,*(BS0LY(K),  K =1 , P ) • , 10 18 ,2 ( /l 8X ,  1018)) 

00005290 

0413 

6146 

I NE  T* 1 

00005292 

0414 

IN IS  = ISTEP 

00005294 

0415 

I NSET=l 

00005296 

0416 

ISTEP*2 

00005298 

0417 

GO  TO  6150 

00005300 

0418 

6147 

I  STEP- INI S 

00005302 

0419 

INSET  =0 

00005304 

0420 

6150 

GO  TO  272 

00005306 

C 

00005308 

0421 

6200 

IF  (NSX.EQ.O)  GO  TO  8100 

00005310 

0422 

6220 

IF  (  lABSCSTX(NSX) ) #GT  • 1000000 1  GO  TO  6500 

00005320 

0423 

6250 

BRO=$TX ( NSX ) 

00005330 

0424 

6270 

STXINSX  )=-BRO-lOCOOOO 

00005340 

0425 

IF  (ISTEP. EO. 0)  GO  TO  6308 

00005390 

0420 

PRINT  6305,  8R0 

00005400 

0427 

6305 

FORMAT!  *0 • ,  *  BR  0  ",110) 

00005410 

0428 

6308 

IF  (ISTEP. NE. 2)  GO  TO  6330 

0000  5*20 

0429 

PRINT  6088,  (STX(INS),  INS=I,NSXl 

00005430 

C 

00005490 

C 

CAPACITY  RULE 

00003500 

94 
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0430 

6330 

NCDsNOOU 

00005510 

0431 

6410 

IF  (ET.EQ.0.0)  GO  TO  6450 

00005512 

l 

0432 

IF  UNSET. EQ.l)  GO  TO  6445 

00005516 

0433 

IF  (INET.EC.I)  GO  TQ  6450 

00005518 

0434 

call  TtMETUNTI 

00005520 

04  35 

ELTN»(INT-lT0)*26.04E-6 

00005523 

0436 

IF  (ELTN.LT.LT)  GO  TO  6450 

00005526 

0437 

6420 

PRINT  6125 «  NOD «  ELTN«  BUB,  BUBS,  1BN00 

00005528 

o 

0438 

IBUB*BUB 

00005530 

0439 

IF  (IBUB.E0.9V99999)  GO  TO  6442 

00005532 

0440 

6430 

PRINT  6135,  (BSOLXt J),J=1,N) 

00005533 

© 

0441 

6440 

PRINT  6145,  (BSOLY(K),  K*lfP) 

00005536 

0442 

6442 

INET-1 

00005538 

0443 

IN  1 S* I  STEP 

00005540 

0 

04  44 

INSET  =1 

00005542 

0445 

7STEP*2 

00005544 

0446 

GO  TO  6450 

00005546 

o 

0447 

6445 

I STEP=INIS 

00005548 

0443 

INSET  =0 

00005550 

0449 

6450 

GO  TO  210 

00005552 

c 

0450 

6500 

IF  (  STX(NSX) .GT. 1000000)  GO  TO  6520 

00005555 

0451 

LX=-STX(NSX)-1000000 

00005560 

0452 

I X  =LX/100 

00005570 

0453 

JX*LX-IX*100 

00005580 

0454 

CX(IX, JX>=0 

00005590 

0455 

KTZ< JXI*KT2< JX)-1 

00005595 

r 

04  56 

GO  TO  6550 

00005600 

0457 

6520 

LX*  S TX ( NS X )  — 1000000 

00005610 

0*»58 

1 X*LX/1Q0 

00005620 

> , . 

0459 

JX=LX-IX*100 

00005630 

046 0 

CX(IX,JX)*0 

00005640 

0461 

F.X( JX)=0 

00005650 

0462 

LQ1*LQ1-1 

00005660 

0463 

6550 

NSX=NSX-1 

00005690 

0464 

GO  TO  6200 

00005700 

. ' 

C 

00005730 

0465 

8100 

1BUB=BUB 

00005740 

0466 

CALL  TIMETUTl) 

00005750 

o 

0467 

ELTl=( ITl-IT0)*26.04E-6 

00005760 

0468 

PRINT  8105,  ELT1 

00005770 

0469 

8105 

FORMAT  <'0*,///lX,  '  ELAPSEO  TIME  IN  SECONDS-*,  F15.8) 

00005780 

0 

0470 

PRINT  8120,  NOO 

00005790 

0471 

8120 

FORMAT  (  • 0# , • TOTAL  NUMBER  OF  NOOES  EXPLORED  «*,I3» 

00005800 

0472 

IF  ( IBUB.6Q. 9999999)  GO  TO  8350 

00005810 

c 

0473 

8130 

PRINT  8150 

00005820 

0474 

8150 

FORMAT  1*0*,  ’NOTE  :  l.  FOLLOWING  X(1,J>  VARIABLES  SHOW  DESIGN • , 

00005830 

l  •  I  TO  WHICH  ACTIVITY  J  IS  ASSIGNED  FOR  J*1  TO  N. •, 

00005840 

9 

2  /7X ,  *2.  IF  EPSILON  EPS  WAS  ASSIGNEO  A  POSITIVE*, 

00005850 

3  *  (NON-ZERO)  VALUE,  THE  SOLUTION  MAY  BE  SUBOPTIMAL. * ,/) 

0000586C 
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